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Abstract 

We are looking for a holographic explanation of nuclear forces, especially the attractive 
forces. Recently, the repulsive hard core of a nucleon-nucleon potential was obtained in the 
Sakai-Sugimoto model, and we show that a generalized version of that model — with an 
asymmetric configuration of the flavor D8 branes — also has an attractive potential. While 
the repulsive potential stems from the Chern-Simons interactions of the U{2) flavor gauge 
fields in 5D, the attractive potential is due to a coupling of the gauge fields to a scalar field 
describing fluctuations of the flavor branes' geometry. At intermediate distances r between 
baryons — smaller than i?KK = O ( 1 ) /M^^ meson but larger than the radius p ~ -RRK/VA'tHooft 
of the instanton at the core of a baryon — both the attractive and the repulsive potentials 
behave as 1/r^, but the attractive potential is weaker: Depending on the geometry of the 
flavor D8 branes, the ratio Ca/r = —V^^^^{r)/V'^'^'^{r) ranges from to |. The 5D scalar fields 
also affect the isovector tensor and spin-spin forces, and the overall effect is similar to the 
isoscalar central forces, V{r) — )■ (1 — Ca/r) x V{r). 

At longer ranges r > -Rkk, we find that the attractive potential decays faster than the 
repulsive potential, so the net potential is always repulsive. This unrealistic behavior may be 
peculiar to the Sakai-Sugimoto-like models, or it could be a general problem of the Nc ^ oo 
limit inherent in holography. 
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1 Introduction 

In recent years, holography or gauge-gravity duality gave us a new approach to hadronic 
physics (see P] for a review). It has been spectacularly successful at explaining many features 
of the quark-gluon plasma such as its low viscosity [2], and there are some interesting results 
concerning the high-density nuclear matter [3l H]. Motivated by this success, the authors 
wanted to apply holography to one of the oldest problems of nuclear physics: The interactions 
between nucleons are very strong, so why isn't the nuclear matter relativistic? Instead, the 
bulk binding energy of the nuclear matter is only 1.7% of Mc^, about 16 MeV per nucleon. 

The usual explanation of this puzzle involves a near-cancellation between the attractive 
and the repulsive nuclear forces: The attractive potential is only a little bit stronger than 
the repulsive potential, and the difference is rather small. For example, in the Walecka's 
mean- field model [S], the attractive potential due to cr- meson field is 400 MeV while the 
repulsive potential due to cj-meson field is 350 MeV; there is also the Fermi motion energy of 
about 35 MeV/nucleon, so the net binding energy is only 16 MeV/nucleon. There have been 
many similar (but more elaborate) models since Walecka, but they all beg the same question: 
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Why is the attractive — interaction between nucleons only a httle bit stronger than the 
repulsive +(\l/7^\l/)^ interaction? Is this a coincidence depending on quarks having precisely 
3 colors and the right masses for the u, d, and s flavors? Or is this a more robust feature 
of QCD that would persist for different Nc and any quark masses (as long as two flavors are 
light enough)? 

The most direct way of applying holography to these issues would be to build a holo- 
graphic model of the bulk nuclear matter. Unfortunately, this approach is troubled by the 
large Nc limit which is inherent in the Holographic QCD. Indeed, even taking the lead- 
ing l/Nf. corrections into account is very hard in HQCD because it requires doing string 
loop calculations on the "gravity" side of the gauge-gravity duality. But for large N^, the 
low-temperature low-pressure phase of the bulk nuclear matter becomes a crystalline solid 
instead of the Fermi liquid for the real-life Nc = 3, and its other properties — such as density 
or the binding energy — could also be quite different. 

Paradoxically, direct holographic modeling works for exotic phases of nuclear matter — 
such as the quark-gluon plasma and maybe the high-density solid phase, if it exist^ — but 
not for the good old nuclei themselves. However, we may still use holography to obtain 
phase-independent features of nucleons and nuclear forces, but relating those features to the 
experimental properties of real nuclei has to be done by some other methods. So in this 
article, instead of trying to model whole nuclei, we focus on a rather humble problem of 
obtaining an attractive two-body nuclear force from the Holographic QCD. 

The first holographic model of a baryon appeared in [S] and [0] in the AdSs x context: 
A D5 brane wrapping the had Nc strings attached to it, while the opposite ends of those 
strings connected to the Nc external quarks at the boundary of the AdSs space. Similar 
"external" baryons were constructed in confining backgrounds in [10]. To make a baryon 
out of dynamical rather than external quarks one needs to add Nf flavor branes to the 
holographic model; usually one takes Nf ^ Nc so the flavor branes act as probes of the 
background created by the color branes. A prototypical model of this kind was constructed 
by Sakai and Sugimoto [TT]: Starting with the Witten's model [T2] of Nc D4-branes on a 



^Similar to the liquid helium solidifying under pressure, the Nc — 3 nuclear matter may also have a 
crystalline high-pressure phase. Although at very high pressures and densities, the nucleons are believed 
to merge into a quark liquid, and it is not clear if the nucleons form a lattice before merging, of if there is 
a direct transition from the nuclear liquid to the quark liquid. If the solid nuclear phase exists at all and 
behaves like a semi-classical crystal, then its structure should not depend much on the Nc and it could be 
modeled using large-iVc methods such as the skyrmion lattices of 6^ , or the holographic instanton lattices of 
[7]. Alas, judging the phenomenological success of such models is rather difficult because the high-pressure 
nuclear matter is hard to study experimentally; the best data comes from modeling neutron stars, and we 
still do not know for sure if their interiors are solid or liquid. 
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circle (with antiperiodic boundary conditions for the fermions to break the = 4 SUSY to 
Af = 0*), they have added Nj D8 and Nf D8 branes. On the gravity side of the duahty, 
the lOD geometry is warped R^''^ x x a cigar, while the D8 and D8 branes connect to 
each other and span R^'^ x x a U-shaped line on the cigar (see figure ([T]) on page [2U]) . 
A holographic baryon comprises a D4 brane wrapping a compact 5*^ and Nc open strings 
connecting this D4 to the flavor D8 branes. To minimize the baryon's energy, the D4 brane 
acting as a baryonic vertex becomes embedded in the D8 branes [T3j, and for Nf > 1 it 
dissolves into an instanton of the U{Nf) gauge theory on the flavor branes. 

Sakai, Sugimoto et al wrote several papers [TU [1^] about properties of such holographic 
baryons, and eventually [16] worked out a repulsive force between two such baryons. But 
they could not get an attractive force because of the accidental Z2 symmetry of the antipodal 
configuration of the flavor branes. To see the connection, note that in the large Nc limit, the 
nuclear forces are dominated by the single-meson-exchange diagrams; the repulsive central 
forces come from exchanges of the vector mesons while the attractive central forces come from 
the scalar mesons. In holography, the 4D vector mesons are modes of the gauge fields living 
on the flavor branes, while the 4D scalar fields are modes of the scalar fields parametrizing 
transverse motion of those branes. In the original version of the Sakai-Sugimoto model, the 
D8 and D8 branes cross the color D4 branes at antipodal points of the 5*^ circle, hence the 
name "antipodal model" . On the gravity side of the gauge-gravity duality, the becomes 
the circular dimension of the cigar, and the combined D8 + D8 branes stretch along the cigar's 
diameter. The two sides of this diameter are symmetric, and this leads to the $ — ?► — $ 
symmetry of the 9D scalar fields parametrizing the transverse motion of the flavor branes; in 
4D terms, this Z2 symmetry flips the signs of all the scalar meson fields, 0j(x) — > —(f)i{x). But 
this symmetry does not affect the gauge fields on the flavor branes and hence the 4D vector 
mesons or the holographic baryons. Consequently, in the antipodal model, the baryons have 
Yukawa couplings to the vector mesons but not to the scalar mesons, and that's why there 
are no attractive nuclear forces but only the repulsive forces. 

In this article, we investigate nuclear forces in the non-antipodal version [17] of the Sakai- 
Sugimoto model. Without the accidental Z2 symmetry, the baryons should have Yukawa 
couplings to both vector and scalar mesons, and indeed we find both repulsive and attractive 
forces. Unfortunately, the attractive forces are too weak and the net force is repulsive at all 
distances, so our model of HQCD is not too realistic. Specifically, at intermediate distances 
r between two nucleons — shorter than the ranges of the 4D Yukawa forces but longer than 
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the size of a baryon's cor^ — both the repulsive and the attractive potentials behave like 5D 
Coulomb potential and scale like But the attractive potential has a smaller coefficient, 

"^attractive 

= ^repulsive = 9 ^ ~ ^ ' '-^'^^ 

where C ^ 1 parametrizes the geometry of the flavor D8 branes: The near-antipodal models 
have ^ 1 and Cajr ^ 1 while the far-from- antipodal models have C ^ 1 ^-^id Cajr ~ |- In 
any case, Cajr < | < 1 and the attractive nuclear potential is weaker than the repulsive. 

At longer distances, nuclear forces are dominated by the 4D Yukawa potentials of the 
lightest mesons with the right quantum numbers, J^'~^ = 1 for the repulsive force and 
J^*^ = for the attractive force, thus 

/ lightest \ / lightest \ 

y repulsive ^ | 6Xp[^— r X TTl^^^^^j. ) ^attractive ^ ^-^Pl~''" ^ "^scalar ) (12) 

In real life, the lightest isoscalar scalar meson cr(600) is lighter than the lightest isoscalar 
vector meson uj{787), so at long distances the attraction wins over the repulsionjfl But in the 
Sakai-Sugimoto models — antipodal or non-antipodal — the lightest scalar meson has more 
than twice the mass of the lightest vector. Consequently, the attractive force has a shorter 
range than the repulsive force, and the net nuclear force is repulsive at all distances. 

We don't know why the meson spectra — and hence the nuclear forces — in the Sakai- 
Sugimoto model are so unrealistic. It could be something peculiar to the model's setup, 
hopefully to be remedied by some future holographic models. But it could also be a general 
problem of the large limit; indeed, the QCD origin of the cr(600) meson is poorly under- 
stood, and it's not clear if for Nc oo it continues to exist or disappears from the spectrum. 
The best way to resolve this issue would be to find the a resonance and its mass in a lattice 
QCD calculation for several values of N^, then extrapolate to Nc — )■ oo. Alternatively, once 
we have several different holographic models, we can compare their predictions for the meson 
spectra in general and for the lightest scalar meson in particular. Either way, this issue will 
have to wait for future research. 



^In all versions of the Sakai-Sugimoto model, the 5D instanton at the core of a baryon has a very small 
size p ^ Rkk/V^ where Rkk is the Kaluza-Klein scale of extra dimensions and A — NcQym S> 1 is the 
't Hooft coupling. On the other hand, the 4D mesons have masses Mmcson ~ ^/Rkk so the Yukawa forces 
have ranges -Rvukawa ~ Rkk S> p. This is quite different from the real-life mesons and baryons where 

Pbaryon ^ 47? Yukawa ■ 

•^Actually, since in real life A^c = 3 ^ oo, the longest-range component of the attractive force is not the 
single-sigma-meson exchange but rather the double- pion exchange. The Yukawa range of this force is 1 /2m^, 
which is significantly longer than l/m^. 
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The rest of this paper is organized as follows. In the next section (§2) we explain the 
problems with the large limit of nuclear physics. First, we explain why large makes the 
nuclear matter solid rather than liquid. Next, we discuss the Nc oo limit of the nuclear 
forces and what happens to the o"(600) meson. Finally, we bring up the issue of separating 
nucleons from other baryonic species such as A. 

Section 3 is a review of the Sakai-Sugimoto model and its antipodal and non-antipodal 
versions. In particular, we derive the effective 5D Lagrangian for the U (2) flavor gauge fields 
(for simplicity we work with two flavors) , then realize a holographic baryon as a lowest-energy 
YM instanton and calculate its mass and radius. Section 4 explains general properties of 
the holographic nuclear forces in the near, intermediate, and far zones; the three zones are 
illustrated in the diagram (jl]) on page |23l We also summarize the calculation by Hashimoto 
et al [16j of the repulsive force in the intermediate zone. 

Section 5 is the core of this paper, that's where we calculate the attractive and repulsive 
nuclear forces in the intermediate and far zones. In §5.1 we derive the effective 5D theory of 
scalar and vector fields living on the flavor branes. We show that the abelian vector and scalar 
fields give rise to 5D Coulomb forces between SU{2) instantons. In the near and intermediate 
zones, both the repulsive potential due to abelian vector and the attractive potential due 
to the abelian scalar have the same dependence, but the attractive potential has a 
smaller coefficient as in eq. (11. ip . In §5.2 we leverage this result to obtain both isoscalar and 
isovector forces between two spinning nucleons at intermediate distances from each other. 
The isovector spin-spin and tensor forces stem from the small overlap between the SU{2) 
instantons implementing the two nucleons and their interactions with the abelian vector and 
scalar fields. Our analysis follows Hashimoto et al [16j, but taking the scalar fields into 
account reduces the isovector forces by the same overall factor (1 — Ca/r) as the net isoscalar 
repulsive — attractive force. Thus, 

Vr,et{r, Jl, h, Jl, J2) = (1 - Ca/r) X 14et('^, h, Jl, J2)[without scalar fields]. (1.3) 

In §5.3 we consider the attractive forces in the far zone. Since in the Sakai-Sugimoto model 
the lightest scalar meson is heavier than the lightest vector meson, the attractive force decays 
with distance r faster than the repulsive force, so the net isoscalar central force is always 
repulsive. We also consider the long-range isovector tensor force due to pions. Although the 
pions are zero modes of the 5D vector fields and have nothing to do with the 5D scalar fields, 
the pion-nucleon coupling depends on the baryon's radius p which is affected by the scalar- 
mediated forces in the near zone. Consequently, the isovector force due to pion exchange is 
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reduced by the overall factor (1 — Ca/r)- Likewise, all other isovector forces in the far and 
intermediate zones are reduced by the same overall factor. 

Our calculation in §5 are based on Yang-Mills approximation for the effective Lagrangian 
for the flavor gauge fields. In section 6 we investigate the validity of this approximation by 
working with a complete non-abelian DBI+CS Lagrangian for the fields on the flavor branes. 
We show that although the SU{2) gauge fields become strong {27ra'J'MN ~ gMN) near the 
center of a baryon, the self-duality of those fields (in four non-compact space dimensions 
of the D8 branes) leads to cancellation of all the higher-order tr(J^^), etc., terms in the 
expansion of the DBI Lagrangian. Consequently, all our calculation in §5 are valid to the 
leading order in 1/A. The leading post-YM effect in 5 D is a small (0(1/A)) correction to 
the self-duahty condition for the SU{2) gauge fields due to abehan vector and scalar fields. 
To see how this correction affects a stand-alone semiclassical baryon, in §6 we minimize the 
DBI -I- CS action of an 5'0(4)-symmetric instanton-like field configuration with a general 
radial profile and show that the minimum is very close to good old YM instanton of the 
same radius as we had in §5. Calculating the non-abelian DBI action involves computing a 
symmetrized trace; this is done in the Appendix. 

Finally, section 7 summarizes our results and makes suggestions for future research. 



2 Limitations of the Nc ^ oo Limit and Holography 

The large N(. limit is inherent in all holographic QCD methods, and this poses a problem 
for the aspects of nuclear physics that are different between the small Nc and the large Nc 
regimes. In particular, the bulk nuclear matter at zero temperature and pressure (but finite 
density) forms a quantum liquid for small — such as real-life N(. — 3 — but becomes a 
crystalline solid for large Nc- 

To see how this works, consider a condensed matter analogy — some atoms which attract 
to each other at long or medium distances but have repulsive hard cores. Semi-classically, 
at zero temperature and pressure such atoms always form some kind of a crystal; it takes 
strong quantum effects to put the atoms into some other phase such as liquid or super-solid. 
Of particular importance is the kinetic energy of the zero-point quantum motion of atoms 
confined to narrow potential wells, 

K (2 1) 

2Matom(well diameter) 2 ' ^ ' ^ 
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or rather its ratio K/U to the potential binding energy U per atom. According to Newton 
Bernardes [18], this ratio is related to the de Bour parameter of the inter-atomic potential 

as 

where Tc is the radius of the atomic hard core and e is the maximal depth of the potential. For 
small de Bour parameters, the quantum corrections to the semi-classical approximation are 
weak and the crystal remains stable at zero pressure. For larger A^, the quantum corrections 
due to kinetic energy become important, and when A^ exceeds a critical value somewhere 
between 0.2 and 0.3 [12], the crystal melts into a quantum liquid]^ For example, helium 
atoms have A^ = 0.306 and hence K/U ^ 1 while neon atoms have A^ = 0.063 and hence 
K/U ~ 0.05; consequently, at zero temperature and zero pressure helium is a quantum liquid 
while neon is a crystalline solid. 

To see how the K/U ratio of the nuclear matter depends on the number of colors, we 
note that in the large limit, the leading nuclear forces are proportional to Nc. Specifically, 
according to Kaplan and Manohar [20] . 



V{r,hj2,J2,J2\N,) = N,^Ac{r)+ A^, x A5(r)(Iil2)(JiJ2) 

+ AT, X AT(r)(Iil2) [3(nJi)(nJ2) - (J1J2)] (2.3) 
+ 0{l/N,). 

for the same A'c-independent radial profiles As, At of the central, spin-spin, and tensor 
potentials. Classically, such potentials would like to arrange a many-nucleon system in some 
kind of a crystal with A'c-independent nearest-neighbor distance ~ 1 fm, while the binding 
energy of a nucleon in such a crystal would be proportional to the Nc- Indeed, all models 
of nuclear matter based on semi-classical models of nucleons form such crystals, for example 
skyrmion crystals of ref. [6]. In the quantum theory, nucleons in such a lattice have zero- 
point kinetic energies ( 12.11) where the well diameter is independent on while the nucleon's 
mass M oc Nc, hence K oc 1/Nc and 

^ oc ^ = ^ (24) 



^Melting releases the individual atoms form narrow potential wells, which significantly lowers their kinetic 
energies ()2.f |) . It also, moves the atoms away from the minima of the attractive potential, which lowers the 
potential binding energy U . The overall effect on the net K — U energy per atom depends on the K/U ratio: 
For low ratios the potential energy is more important and the crystal is stable, but for high ratios lowering 
the kinetic energy becomes advantageous and the crystal melts. 



7 



We may estimate the coefficient of tliis proportionality using tlie de Bour parameter A^. Tlie 
maximal depth of the central potential between two nucleons is about 100 MeV for A'^^ = 3, 
so we take it to be e ~ A^c x 30 MeV for large Nc- Likewise, we take the nucleon mass to be 
Mtv ~ A'c X 300 MeV and hard-core radius rc ~ 0.7 fm regardless of Nc. Consequently, 



and hence liquid nuclear matter for A"c < 8 and solid nuclear matter for A^c ^ 8. 

The numerical coefficient in eq. (I2.5P and hence our estimate A"^"^'* ~ 8 for the dividing 
line between liquid and solid bulk nuclear matter (at low pressures and temperatures) should 
be taken with a large grain of salt. Also, the transition between liquid nuclear matter for 
A'c = 3 and crystalline nuclear matter for large Nc may go through some exotic phases at 
intermediate values of Nc, perhaps something like a quantum supersolid, perhaps something 
more exotic without known condensed-matter analogues. But regardless of the details of 
this transition, in the large Nc limit the potential energy of interacting near-static nucleons 
becomes much larger than the nucleons' kinetic energies, and the bulk nuclear matter at 
T = 0, P = conditions becomes a conventional semi-classical crystal. The structure of 
such crystals can be modeled holographically — and indeed there is active research in this 
direction (for instance [7J) — but we have no experimental data to compare to the models 
because real-life nuclei with A^ = 3 are liquid rather than solid. 

Meanwhile, instead of trying do build holographic models of complete nuclei we focus on 
holographic models of the nuclear forces. But even at the level of the two-body forces, the 
large A^ limit maybe different from the real-life case of just 3 colors. Of particular concern 
is the isoscalar attractive force due to exchanges of the cr(600) scalar mesons between the 
nucleons. In real life, this is a major component of the net attractive force — especially at the 
medium-long distances between the nucleons — but in the large Nc limit this component may 
weaken or disappear because the a (600) meson itself may become heavier or even disappear 
from the scalar meson spectrum. 

The (7 (600) (also known as /o(600)) is the lightest isoscalar true-scalar meson. In real 
life, it appears as a very broad resonance of two pions — so broad that its central mass is 
somewhat controversial and different experimentalists locate it anywhere between 400 meV 
and 700 MeV, and sometimes even higher, cf. references in the Particle Data Group's listing 
[2T] . But the real controversy about the (j(600) resonance is its physical origin. Unlike the 
heavier = 0"*", J*"^ = O"*""*" mesons /o(980), /o(1370), etc., the a (600) meson does not exist 




2 



K 



45 



(2.5) 



8 



in the non-relativistic quark modeo so for many years R. L. Jaffe and others [221 1231 1211 [25] 
were claiming that the ex (600) is not a true qq meson but a qqqq tetraquark. Specifically, it's 
a molecule-like bound state of two pions which exists because the p-meson exchanges in the 
t-channel induce an attractive s-channel force between the pions. If this claim is true, then 
the a resonance goes away in the large Nc limit because the forces between pions become 
weak as 1/Nc. 

But many other authors (see [26] for a sample) identify the a (600) with the a field of 
the linear sigma model of the chiral symmetry breaking. Or rather, the massive (j{x) field 
parametrizing fluctuations of magnitude of the symmetry-breaking VEV gives rise 

to primordial sigma quanta, while the real sigma mesons cr(600) are quantum mixtures of 
those primordial quanta with the Ivrvr) states (and to lesser extent with the other J"^ = 0^, 
J^^ = 0++ mesons). From this point of view, the non-relativistic quark model is irrelevant 
because the quarks do not become non-relativistic until after the chiral symmetry has already 
been broken. Indeed, the NRQM does not see that the pions are (pseudo) Goldstone bosons, 
so the fact that it does not see the sigma meson at all is simply another limitation of the 
NRQM as far as the chiral symmetry breaking is concerned. If this point of view is right, 
then the sigma meson exits for all Nc- For large A^^ limit, this meson is mostly a quantum 
of the (y{x) field — its mixing with Itttt) and other states becomes weak — and it's a narrow 
resonance rather than a broad hump we have for A^^c = 3, but it remains a dominant resonance 
in the l'^ = 0~^, J^^ = 0^^ rni channel, and its mass should not be too different from the 
real-life 600 MeV. 

The other mesons — scalar or vector, isoscalar or isovector — are unlikely to be disturbed 
by the large limit, so their contributions to the nuclear forces would be similar to real-life 
QCD. IF the ct(600) meson remains in the spectrum in the large A^^^ limit and if its mass 
remains similar to the real-life 600 MeV, then the entire nuclear potential (12. 3p for Nc ^ oo 
would be similar to what it is in real life, except for the overall factor N^. In particular, 
the net central potential Vc{r) would be repulsive at short distances (the hard core) but 

^In the non-relativistic quark model, all 0"'""'" mesons have 5 = 1 and L = 1. Consequently, the lightest 
0++ meson should be heavier than the lightest 1 mesons p(770) or w(787) that have S ~ 1 but L ~ 0. 
Depending on the assumptions one makes about the forces between the quark and the antiquark, this 
argument identifies the lightest true qq meson with = 0+ and J^<^ = 0++ as either /o(980) or /o(1370). 
In any case, the cr(600) resonance is way too light to be a p-wave qq state, so it has to be something else. 
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attractive at medium and long distances: 



K(r) 




iVc — )• oo limit with a light a meson 



(2.6) 



On the other hand, if the (t(600) meson disappears from the spectrum for large N^-, or 
if it becomes heavier than the lightest vector meson, then the dominant attractive force 
would become shorter-ranged than the repulsive force, and the net force at medium and long 
distances would be repulsive rather than attractive: 



In this scenario, at large Nc the nuclear force is repulsive at all distances, and there are no 
bound nuclei at all, liquid or crystalline. 

So what really happens to the sigma-meson and to the nuclear forces at large N^"^ The 
best way to settle this controversy would be to find the a resonance and its mass in a lattice 
QCD calculation for several values of N^.. Such a calculation would require a realistic pion 
mass (unlike most present-day lattice calculations extrapolating from m,r > 350 MeV) and 
rather large lattices to distinguish the sigma resonance from the two-pion continuum, so it 





Nc ^ oo limit without the a meson 



r 



(2.7) 
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may be too hard for the present-day computers. But thanks to the Moore's Law, finding the 
a resonance on a lattice should become possible in a not-too-distant future. 

Alternatively, we may try to resolve the issue using holography. Although a holographic 
model of real QCD — or rather, of QCD with large N^. — is yet to be constructed, several 
known models seem to be qualitatively similar, so we can compare their predictions for the 
meson spectra in general, and for the lightest true scalar meson in particular. However, 
the models that seem qualitatively similar to QCD may not be similar enough, and their 
predictions could be widely off target. Indeed, the predictions of different models have 
turned out to be quite different from each other. For example, in the Sakai-Sugimoto model 
which we use in this article, the lightest true scalar meson is more than twice as heavy as the 
lightest vector meson. Consequently — as we shall see in painful detail in section 5 — the 
net nuclear force is everywhere repulsive and looks like (12. 7p rather than like (12. 6p . On the 
other hand, in the highly-non-antipodal version of the Dymarsky-Kuperstein-Sonnenschein 
model [27], the lightest J^^ = 0++ meson is much lighter than any other mesons (except the 
pions) [29]. However, this lightest scalar is a pseudo-Goldstone boson of the approximate 
conformal symmetry of the flavor sector, so it is not clear how much attractive force it can 
mediate. As of this writing, it is not clear if the net nuclear potential in this model looks 
like the real-life potential (12. 6 p or like the everywhere-repulsive potential (12. 7p we calculate 
in this paper for the Sakai-Sugimoto model. 

But suppose tomorrow somebody discovers a holographic model of the real QCD and — 
miracle of miracles — it has a realistic spectrum of mesons, including the cr(600) resonance, 
and even the realistic Yukawa couplings of those mesons to the baryons. Even for such a 
model, the two-body nuclear forces would not be quite as in the real world because the 
semi-classical holography limits A^^^ — )■ oo, A — )■ oo suppress the multiple meson exchanges 
between baryons. Although in this case, the culprit is not the large number of colors but 
the large 't Hooft coupling A = Ncgyj^. 

Indeed, from the hadronic point of view, nuclear forces arise from the mucleons ex- 
changing one, two, or more mesons, and in real life the double-meson exchanges are just 
as important as the single-meson exchanges. In particular, since the lightest mesonic state 
with J*^ = 0"^, J'"^ = O^"*" quantum numbers is a pair of un-bound pions, the longest-range 
isoscalar attractive force between nucleons comes from exchanges of two pions rather than 
of any single mesons. In holography, the single-meson exchanges happen at the tree level of 
the string theory while the multiple meson exchanges involve string loops (fc — 1 loops for k 
mesons), and the loop amplitudes are suppressed by the powers of 1/A relative to the tree 
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amplitudes. 

Naively, one would expect the loop amplitudes to carry additional factors oil/Nc (which 
is dual to the string coupling) rather than 1/A, or maybe both 1/A and factors, but the 
naive power-of-A^^c counting does not work for loop amplitudes involving baryons made of N^. 
quarks]^ Indeed, in honest QCD with a large number of colors, the multi-meson-exchange 
contributions to the non-relativistic effective potential for the baryons are not suppressed by 
powers of 1/Nc [27]. However, the extra powers of Nc due to quarks in a baryon are not 
accompanied by the extra powers of A, so in holography, the contributions of the multiple 
meson exchanges are suppressed, albeit by powers of 1/A rather than l/N^. 

To see how this works in a general holographic model of QCD with ^ Nf, note that 
such a model starts with a string-theoretic construction where the colors and the flavors live 
on separate branes. For large and large A, the color branes become black branes producing 
curvature and fluxes through the bulk, which provide a non-trivial background for degrees of 
freedom living in the bulk itself as well as on the flavor branes. The bulk degrees of freedom 
are dual to the pure-color sector of QCD (glueballs, etc.), while the vector and scalar fields 
living on the flavor branes are dual to the qq mesons. The flavor fields have rather weak 
couplings to each other: in 5D terms. 



^5D,flavor - (2.8) 

SO the 4D mesons — which are modes of the 5D vector and scalar fields with wave functions 

— 1/2 

il) ~ Rkk — have couplings to each other of the order 

/ 

/ 

/ 

/ 

/ 

1 

9MMM - (2.9) 



\ 
\ 

\ 

\ 

\ 

\ 



A holographic baryon is made from some brane spanning only the compact dimensions that 
is connected to the flavor branes by N^. strings, although this construction is often equivalent 
to an instanton of the 5D flavor gauge fields. Consequently, the baryon-meson coupling is 



^The authors thank Alexei Cherman and Thomas Cohen for bringing this fact to our attention after we 
made a mistake in the first version of this paper. 
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enhanced by an extra factor of N^, 




9mbb ~ X 



1 %/iv: 



(2.10) 



At the tree level of the baryon-meson theory, scattering of two baryons proceeds through 
a single-meson exchange, which produces a 0{Nc/X) amplitude. 





J tree 2 



A 



(2.11) 



At the one-loop level, there are two types of diagrams, the triangle diagrams such as 





^ 9mbb ^ 9mmm ~ -jj (2-^2 




and the box and crossed-box diagrams 




+ 





\ / 
\ / 
\ ✓ 

X 

/ \ 
/ \ 
/ \ 





(2.13) 
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with amplitudes 

^° - 9mbb - -f ■ (2.14) 

that carry an extra power of N^. However, Banerjee et al showed [27j that for non-relativistic 
baryons, the box and the crossed-box diagrams almost cancel each other from the effective 
potential between the baryons, with the un-canceled part having a lower power of the N^.- 
Banerjee et al did not pay any attention to the powers of A, but clearly the un-canceled 
sub-leading terms in the box and crossed-box diagrams cannot carry higher powers of the 
't Hooft coupling than the leading terms fl2.14p . thus 

N 1 

^uncanceled ~ ~ ~ -A^'^". (2.15) 

In other words, the contribution of the double-meson exchange and other one-loop processes 
to the 2-body nuclear potential carries the same power of but is suppressed by a factor 
1/A compared to the tree-level singe-meson exchange. 

To be precise, the large A limit suppresses exchanges of the un- bound meson pairs but 
not of the meson-meson resonances — which become narrow (because of weak q^mm) ^^"^ 
act as single mesons exchanged between the two baryons. In particular, in the isoscalar O^"*" 
channel that gives rise to the dominant attractive force between nucleons, the A — oo limit 
suppresses the contribution of the unbound two-pion continuum, but it replaces it with a 
discrete set of /o resonances. In a good holographic model of QCD (which alas has not been 
found yet), the overall strength of the O^"*" channel should be similar to the real QCD, so it 
would produce a similar isoscalar attractive force at short distances. However, the range of 
this attractive force would be significantly shorter: Instead of decaying with distance like 
exp(— 2m^r) as in real life, the holographic attractive force decays as exp(— mor) where mo 
is the mass of the lightest isoscalar O^"*" meson, presumably cr(600 MeV). 

In principle, the isoscalar 1 channel that gives rise to the dominant repulsive force 
suffers from similar corrections in the A — oo limit. But in practice, the strongest and the 
longest-range contribution to this channel comes from exchanges of a single a; (787) meson, 
so suppressing the multi-meson exchanges in this channel would not make a qualitative 
difference. Thus altogether, the net effect of large 't Hooft coupling on the central nuclear 
potential — besides the overall 1/A factor — is the shortening of the attractive tail at long 
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distances: 



A X V,{r) 



blue: real QCD, A ~ 1 

red: best possibility for holographic QCD, A ^ 1 




(2.16) 



r 



However, this optimistic picture presumes a holographic model of QCD that correctly repro- 
duces (a) the overall strength of the isoscalar 0"'"'*' and 1 channels, and (b) the mass spectra 
of vector and scalar mesons, especially the masses of the lightest O"*""*" and 1 mesons (7(600) 
and (X'(787). But thus far, no known model satisfies these requirements, not even approxi- 
mately, so the nuclear forces they produce could be much more different from the real life 
than fl2.16p . In particular, the nuclear force we calculate in this paper for the Sakai-Sugimoto 
model turns out to be everywhere repulsive: 

K(r)/iV, 




(2.17) 



Now let's go back to the large Nc limit — in holography or in honest QCD — and consider 
yet another general problem with baryons made from many quarks: How to separate the 
nucleons with J = J = ^ from the other kinds of baryons such as A with J = J = |? 
In real life, there is a large mass gap between the nucleons and the A baryons — almost 
300 MeV — but for large this gap shrinks as l/N^. At the same time, the two-baryon 
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potential grows like N^, so for large it becomes stronger than the gap. Consequently, two 
interacting nucleons may "forget" their individual spins and isospins and mix up with other 
baryonic species such as A. In fact, for large Nc there is a whole lot of baryonic species 
with I = J ranging from | (for odd Nc) or (for even Nc) all the way up to Nc/2, and 
a strongly-interacting nucleon might mix up with all of them. While such mixing would 
not affect the isoscalar spin-blind central force between two baryons, it might significantly 
enhance the isovector spin-spin and tensor forces. 

Therefore, comparing the two-baryon forces in the large Nc limit to the real-life two- 
nucleon forces is rather tricky. One has to carefully keep track of the spin and isospin 
degrees of freedom of the two baryons, expand the interaction Hamiltonian into central, spin- 
spin, and tensor forces as in eq. (12. 3p . and then compare the radial profiles Ac{r), As{r), 
and Axir). Moreover, the spin and isospin degrees of freedom require quantum mechanical 
treatment because semi-classically, we do not get definite spins or isospins even for stand- 
alone single baryons. Instead, we get skyrmions, or instantons, or some other kind of solitons 
with a definite orientation of the 5'f/ (2) isospin relative to the S'f/(2)spin; in quantum terms, 
they become superpositions of baryons with all possible J = J= |,|,...,oo. Consequently, 
a force between two such semiclassical baryons is not a force between two nucleons but rather 
a superposition of forces between different baryonic species. 

This problem affected the first holographic calculation of the nuclear forces by K. Y. Kim 
and I. Zahed [30]. Their baryons were semiclassical instantons in the Sakai-Sugimoto model, 
so instead of definite \I, Iz, J, Jz) they had a definite direction n in S*^ = 5'f/(2)isospin x 
S'f/(2)spin/5'[/(2)common- Consequently, Kim and Zahed[30] found that the force between two 
baryons depends on the angle between rii and 112 — it was attractive for some angles and 
repulsive for other — but they could not interpret this angular dependence in terms of the 
isovector spin-spin and tensor forces. By comparison, Hashimoto, Sakai, and Sugimoto [16j 
made a similar calculation using properly quantized collective coordinates for each instanton. 
Consequently, they obtained the force between two nucleons rather than some mixed-up 
baryons, and they could see how this force depends on each nucleon's Iz and Jz- In particular, 
they saw that at medium-short distances, the net force between two nucleons is always 
repulsive. Evidently, the attraction Kim and Zahed saw for some relative orientations of 
semiclassical baryons happens only for high spins and isospins, but not for nucleons with 
I = J = i 

On the other hand, the analysis of Hashimoto et al was limited to the first-order pertur- 
bation theory for nucleons that are far enough from each other to avoid the strong mixing 
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of spins and isospins. This approach will not work for the hard-core region at very short 
distances where the interactions are much stronger than the gaps between states of the 
individual baryons. In the hard core, the semiclassical analysis of Kim and Zahed might 
work better than the perturbative expansion of Hashimoto et al, although comparing the 
semi-classical large- A^^c results to the real-life nuclear forces might be problematic. 

To summarize, the large Nc limit of nuclear physics suffers from three major problems. 
The third problem of baryon mixing is only technical, and it can be solved — at least 
for the medium and long distances between the nucleons — by following Hashimoto et al 
rather than Kim and Zahed. But there are no ways around the first problem of different 
phase structures of nuclear matter with A''^ = 3 and with ^ oo. Even at high pressures 
and densities, there is a difference: For A^^^ = 3, squeezing nucleons together makes them 
merge into a quark liquid, while for Nc ^ oo the nucleons always retain their individual 
identities and a would-be quark liquid suffers from the "chiral density wave" instability 
It is possible that at some intermediate pressures and densities the A'c = 3 nucleons form a 
crystal — just like helium solidifies at high pressures — before merging into a quark liquid. 
If this intermediate- pressure phase of real nuclear matter is ever observed in a lab, or can be 
reliably shown to exist in some exotic but observable places like inferiors of neutron stars, it 
would be very interesting to compare its properties to the holographic models. Until then, 
we can only speculate. 

Finally, the second problem — concerning the fate of the a (600) resonance in the large 
Nc limit and its effect on the attractive nuclear force — is solvable in principle, but it has 
not been solved yet. In holography, this problem is aggravated by using QCD-like models in 
lieu of the presently unknown holographic dual of the real QCD. The meson spectra of such 
models are not quite realistic; for example, in the Sakai-Sugimoto model (both antipodal 
and non-antipodal versions) there is no a resonance and the lightest scalar meson has more 
than twice the mass of the lightest vector meson. Consequently, we shall see in section 5 
that in this model, the attractive force is both weaker and shorter-ranged than the repulsive 
force, so the net nuclear force is always repulsive. This could be a peculiar failing of the 
Sakai-Sugimoto model, or it could be the general problem of holography or even of the large 
Nc limit. Hopefully, future research will resolve this issue. 
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3 Baryons in the Non-antipodal Sakai-Sugimoto Model 



The Sakai-Sugimoto [TT] model is based on placing a set of Nj D8 and anti-D8 flavor branes 
into the gravitational background of A^^^ coincident near-extremal DA branes [12]. We take 
the Nc ^ Nf limit, so the flavor branes are treated as probes. The color D4 branes span the 
Minkowski spacetime an a compact circle of radius R in the direction X4. In the holographic 
limit, these branes become a background comprised of the following metric, RR 4-flux, and 
dilaton: 



u 



R 



DA 



3/2 



-dt^ + Sijdx'^dx^ + f(u)dxl 
3/4 



R 



DA 



U 



3/2 



du'' 

L7R 



+ u'^dnl 



X €4 



where 



R^ 



DA 



9s 



u 



R 



DA 



(3.1) 



f{n) 



1 



Ma 
u 



(3.2) 

'a' is the string 



V4 is the volume of the unit sphere Vt^ and €4 is its 4-volume form, 1^ - 
length, and Qs is the string coupling. 

The manifold spanned by the coordinates u and X4 has the topology of a cigar with tip at 
Wmin = ^^a; the flavor branes span a continuous line on that cigar, see figure ([1]). In order to 
avoid a conical singularity at the tip of the cigar, the radius R of the X4 circle has to satisfy 



47r 

2ttR = — 
3 



R^ 



DA 



UK 



1/2 



(3.3) 



Consequently, the Kaluza-Klein scale of the model is 



1 Q 1/2 

^DA 



while the confining string tension [31] is 



(3.4) 



- str 



27r^2 



u= 



1 



2nP \R 



Ma 



i-DA 



3/2 



(3.5) 



Besides the line on the (m, X4) cigar, the D8 branes span the Minkowski spacetime and 
the compact S'^ sphere. All Nf branes are coincident, and the action of this brane stack has 
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a Dirac-Born-Infeld (DBI) term and a Chern-Simons (CS) term, 



'S'ds — 'S'dbi + '5'cs- (3.6) 

The DBI term is 

•S'dbi = Tg j (fx e'^ str (^^/- det{gmn + 2-Ka'J^„in)^ (3.7) 

where Tg = (27r)^®/~^ is the D8-brane's tension, Qmn is the nine-dimensional induced metric 
on the branes. 



u 



3 ,.3 , .2' 

M\^„,x I R\'^ 1 I du 



f{n) + 



3 



dx\ + ( ~ ) u^dVtl, 



(3.8) 

J-'mn is the U{Nf) gauge field strength on the worldvolume of the D8-branes, and str is 
the symmetrized trace over the fiavor indices. The 9D Chern-Simons term involves the 5D 
gauge fields and the RR fiux on the 5"^; after integration over the 5"^ it becomes the usual 
5D CS term. 

In the geometry of the cigar, the D8 and the D8 branes cannot continue indefinitely 
in the negative u directions, and they have no place to end. Instead, the branes connect 
to the antibranes and form continuous U-shaped lines, whose geometry satisfies classical 
equation of motion for the DBI action fl3.7p . The solutions form a family parametrized 
by the uq > ma, the lowest point on the brane stack. Figure ([1]) below illustrates two such 
solutions: The physical meaning of the uo parameter is the "string endpoint mass" of a quark 
defined such that the mass of a stringy meson involving a non-spinning string of length L is 
^meson = ^str-^ + 2m^^^. In terms of the brane geometry, this endpoint mass is 



1 n . . 1 



no 



mr = ^ I V^.iu = — 1^ f-^'Hu) in. (3.10) 

The original Sakai-Sugimoto model assumed the antipodal configuration of the fiavor branes, 
as shown on the right side of the figure ([T]); for this configuration Uq = U\ and m^^^ = 0. But 
comparing holographic predictions to for the meson and baryon masses to the experimental 
data shows that non-zero end-point mass gives a better fit. Likewise, the width of a meson 
due to its decay into two mesons [33], the best fit to Regge trajectories, and matching baryonic 
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X 



X 



Figure 1: The figure on the right is the generahzed non-antipodal configuration. The figure 
on the right describes the hmiting antipodal case L = irR, where the branes connect at 

Uo = Ua. 

properties [13] , all favor non-zero end-point masses and hence non-antipodal models depicted 
on the left side of the figure ([1]). 

In section 6 we shall consider the full DBI -|- CS action for the fiavor branes, but most of 
our analysis will be based on the weak-gauge-field approximation 

^DBib,^] = S'^M + Sym + 0{T^) (3.11) 

where 

is the action for the metric with switched-off gauge fields. To simplify this action — as well 
as the YM action for the gauge fields — it's convenient to change the cigar coordinate along 
the branes from u or X4 to z defined according to [13] 

u = UA{e + Cz'f'\ C = -. (3.13) 

The reason for this choice is that both Xi{z) and u{z) are single-valued smooth functions 
(unhke the x^^iu) function which is double- valued and singular aiu — Uo). Unlike in [TT], here 
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z is dimensionless; its values run from — oo to +00. The parameter is 1 for the antipodal 
model while the non- antipodal models have C > 1- 

The Yang-Mills part of the action (13. lip is obtained by expanding fl3.7p to the lowest 
non-trivial order in gauge fields. In terms of the z coordinate, 

^YM = ^d'xdz tr(iMz; Q^l^ + Mlk{z- Q^l^j (3.14) 

where 

«^ = ^^'^(^yL^' ^ = ^-al-M. Kz-A) = (l + 2;2)-V3, k{z-l) = l + z\ (3.15) 

and the full expressions for h{z\ Q and k{z] Q with C > 1 are given in equation (4.2) of [13]. 
For the case of two flavors, the U (2) gauge fields can be decomposed as 

= ^?U(2) + \^U{1) = M'"" X + i^"" X 12X2 (3.16) 

where A™ are the SU{2) gauge fields, are the Pauli matrices, and A"^ is the abelian gauge 
field in U{1) C U{2). For a baryon, the SU{2) gauge fields form an instanton in the 4 space 
dimensions (1,2,3,2;), while Chern-Simons coupling induces the abelian electric field Aq. 
After rescaling the coordinates and truncating equations of motion to the leading terms in 
1/A, the gauge fields take form [161 [13] 



Am{x\z) = iv{^) X gdMg ^ 

e 



g{x,z) = (« = 1,2,3), (3.17) 



AC ie + ' 

Acl _ 4CI _ r, 
^0 ~ — ^1 

where M, = 1, 2, 3, z and z = g^lzs x j^- Note that there is a critical value for (, namely, 
C = (5/8)^/'^, but this critical value is unphysical — all the Sakai-Sugimoto models have 
C > 1 since the branes cannot go below the tip of the cigar. Substituting the fields (I3.17P 
into the YM + CS action and minimizing with respect to the baryon radius p gives the 
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baryon's classical mass 



Note that similar to the antipodal case, the baryon's radius scales as A~^/^, only the numerical 
coefficient is different for ( ^ 1. 

By repeating the analysis of pL4j for the non-antipodal models, Seki and Sonnenschein [13] 
obtained the mass spectrum for baryons with different I = J and radial quantum numbers, 
the mean-square charge radii, the isoscalar and isovector magnetic moments, etc., etc., as 
functions of ( and Ma. The found the best fit to experimental data obtain for an un-physical 
( ^ 0.942, which indicates that the generalized Sakai-Sugimoto model is not a very accurate 
description of real-life baryons. Nevertheless, in this article we shall stick with the non- 
antipodal Sakai-Sugimoto models simply because its the only model of holographic QCD we 
know in detail. 



4 Summary of the Repulsive Force 

In real life, the nucleon has a fairly large radius compared to the ranges of mesonic Yukawa 
forces (except pion's), -Rnucieon ~ 4/Mpineson- But in the holographic nuclear physics with 
A ^ 1, we have the opposite situation: While the meson masses are 0(Ma) and the 4D 
Yukawa forces have 0{1/Ma) ranges, the baryon has a much smaller radius -Rbaryon ~ 
A^^/^/Ma, cf. eq. f l3.18p . Thanks to this hierarchy, the nuclear forces between two baryons 
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at distance r from each other fall into 3 distinct zones: 




(4.1) 



In the near zone r < -Rbaryon ^ (V^a); ^^e two baryons overlap and cannot be approx- 
imated as two separate instantons of the SU{2) gauge field; instead, we need the ADHM 
solution of instanton number = 2 in all its complicated glory. On the other hand, in the 
near zone, the nuclear force is five- dimensional: the curvature of the fifth dimension z does 
not matter at short distances, so we may treat the U{2) gauge fields as living in a flat 
5D spacetime. To leading order in 1/A, the SU{2) fields are given by the ADHM solu- 
tion, while the abelian Ao{x,z) is the 5D Coulomb field coupled to the instanton density 
(l/327r^)e°^'^*^^ tr(FR-iFjv/Ar)ADHM- Unfortunately, for two overlapping baryons this density 
has a rather complicated profile, which makes calculating the near- zone nuclear force rather 
difficult. 

The far zone r > (1/Ma) ^ -Rbaryon poses the opposite problem: The curvature of the 5D 
space and the z-dependence of the gauge coupling becomes very important at large distances. 
At the same time, the two baryons become well-separated instantons which may be treated 
as point sources of the 5D abelian field A^. In 4D terms, the baryons act as point sources for 
all the massive vector mesons A!^{x) comprising the massless 5D vector field A^{x, z), hence 



23 



the nuclear force in the far zone is the sum of 4D Yukawa forces, 

/\r2 p-mnr 

Vir) = ^El^«(- = 0)l'x^ (4-2) 

n 

where m„ = 0(Ma) are the vector meson's masses and ipniz) are their wave functions in the 
curved fifth dimension. At the inner edge of the far zone, all the 4D vector mesons contribute 
to the potential (14. 2 p but for larger distances, the lightest vector meson becomes dominant. 

In the intermediate zone -Rbaryon r (1/Ma), we have the best of both situations: The 
baryons do not overlap much and the fifth dimension is approximately flat. At first blush, 
the nuclear force in this zone is simply the 5D Coulomb force between two point sources. 

This behavior of the repulsive potential suggests that the intermediate zone of the 
holographic nuclear force corresponds to the repulsive hard core of the real-life nucleons. 

The real-life hard-core repulsion has both isoscalar and isovector components of compa- 
rable strengths, but the potential (14. 3 p is purely isoscalar. The reason for this discrepancy 
is that the point-source approximation of holographic baryons is too crude for the inter- 
mediate zone where two instantons of size p = -Rbaryon have 0((i?baryon/''^)^) effects on each 
other. Moreover, since the size of a stand-alone baryon is a compromise between two sub- 
leading effects — the U{1) Coulomb repulsion and the z-dependence of the gauge couplings 
Kh{z) and Kk{z) — the baryons are linearly sensitive to anything affecting the leading SU{2) 
gauge fields. Hence, the overlap between two baryons gives rise to an additional nuclear 
force of strength comparable to (14. 3p . and since the overlap depends on the baryon's relative 
isospins, this extra force has an isovector component. 

To properly account for the baryon-baryon overlap, Hashimoto, Sakai, and Sugimoto [T6] 
(and also Kim and Zahed [30]) set the SU{2) gauge fields to the self-dual ADHM solution of 
instanton number = 2. The instanton density I{x,z) = (l/327r^)e°'^^^^^tr(Fi^j^FA,fAr)ADHM 
of this solution deviates by 0((i?baryon/'")^) from the sum of two separate instantons, and 
that has two effects: (A) the U{1) Coulomb energy is significantly different from (14. 3p . and 
(B) the width of the instanton density in z direction is different, which changes the SU{2) 
field's energy since the gauge coupling depends on z. Somehow, the two effects cancel out 
from the isoscalar components of the hard-core potential, but they do give rise to isovector 
forces of comparable magnitude. Specifically, for baryons of spin = isospin = | — i.e., for 
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two nucleons — Hashimoto et al obtained 

V{r,n) = ^xlx(l_ + y(Ii-l2)(n-Ji)(n.J,)) = Q (4.4) 

for the intermediate-zone distances r. Note that the isoscalar component of this hard-core 
potential is precisely as in the naive eq. (14. 3p . it's the isovector component that has really 
needed all the hard work. 



5 Attractive Forces in the Non- Antipodal Model 
5.1 5D Scalars and their Interactions. 

In 4D, the attractive forces between two baryons emerge from exchanges of virtual mesons 
with even spins and positive parity, especially the true scalars 0"*". In the holographic theory, 
the baryons are instantons of the 5D non-abelian gauge fields, while the 4D scalar mesons 
are modes of the 5D scalar field $(x), so to get an attractive nuclear force we need a 5D 
scalar- vector coupling of the form 



^5D D Jd'x^xtT{Flj^). (5.1) 



In the Sakai-Sugimoto model, the 5D scalar $(x) describes deviations of the D8 brane stack 
from its equilibrium position in the {u,x^) plane/l For the non-antipodal version of the 
model, such deviations have a 6u component, and since the local 5D gauge coupling depends 
on u, we get (5m($) x tr(F|j^) interactions as in eq. (15. ip and hence the attractive force. 
Unfortunately, for the antipodal model worked out by Sakai and Sugimoto themselves, the 
equilibrium brane stack lies along the radius = const (or rather two opposite radii), so the 
first-order deviations are in the x^ direction only and have no 6u component. Consequently, 
these is no linear $ x tr(F|j^) coupling in 5D — in fact, there is exact Z2 symmetry $ — )■ — $ 
which forbids it — and that's why there is no attractive nuclear force in the antipodal model. 

In this section, we shall derive the effective 5D Lagrangian — including the crucial 
coupling (15.11) — for the non-antipodal model, and then use it to derive the attractive 
nuclear force for the intermediate distance. Our first step is a precise definition of the scalar 



'"In addition to the isosinglet scalar $(a;) which describes the motion of the whole D-brane stack, there are 
also isotriplet scalar fields ^"(x) which describe the relative motion of the two D8 branes. For the moment, 
let us focus on the isosinglet $, we shall return to the isotriplets later in this section. 
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field in terms of tlie D8-brane stack deviation from its equilibrium position. Using 

some kind of a non-singular coordinate w along the brane stack, we define: 



In equilibrium : u = u{w), = x'^{w) (5.2) 

deviation: u{w,x^) = u{w) + 7ra'$(w,x^), (5.3) 
x^{w,x'') = x^{w) - l3{w) X TTa'^{w,x''), (5.4) 

where fi(u.) = ^"ff ^^g^'' ■ (5.5) 

g44{u) X [dx^/dw) 

The last formula here assures that to first order in $, the deviation of the stack is locally 
perpendicular to the stack itself. 

In eqs. fl5.2H5.5p w is a generic non-singular coordinate along the un-perturbed brane 
stack, for example z = [{u^ — uD/uouj^]^^"^ (but not the original z which would be affected 
by the deviation field $). However, to simplify the 5D notations we would like to have the 
same metric for all five dimension, g^w = gii, at least for $ = 0. This calls for 

dw'' = f{u){dx\w)f + -^±-(du{w))\ (5.6) 

which together with the brane equilibrium equation 



du J u^Piu) f {u) - ulf {uq) 



gives us 



du _ u^f{u) - ulf{uo) dx'^ _ ul^J fiyUo) 



dw ]l u^Rf)^ ' dw u^f{u) ' 

implicitly defining the w coordinate. We could not solve these equations analytically, but 
fortunately we would not need the explicit formulae in this paper. All we will need to know 
is that 

(0^3 _ c \ 
1 + ^g^^ X (Maw) + 0{{Mj^w)^)j. (5.9) 
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For the above definitions, the 5D metric for = {x^,w) becomes 



,2 I'u + Tia'^Y^^ 
d4d = I ^ I X 



R 



Mg/(Mo) 

- 8^ X dw' + 0($^) 
u 



(5.10) 



while the radius and the dilaton depend on $ according to 

radius[5^] = R'jiXu + na'^y^\ = ( ^+^^'^ \^'\ (5.11) 

Consequently, expanding the DBl action 

^DBi = T8y"d^xe-'^VoI(^^)str(det((?^,^ + 27ra'J-Miv))'^' (5.12) 
to the second power in J^nm and Om^ and to the first power in $ itself, we get 

5'dbi = d'^xdwl const + £kin + '^mt), (5.13) 



216773 

(( using the r^*^^ to contract the 5D indices, )) 
Ant = ^(-3$ X ItriTp + 5$ X tr(J-^^)) + 0{<^'). (5.15) 

Thus far, we have focused only on the isosinglet scalar field $ describing the common 
motion of the two flavor D8 branes but ignored the isotriplet scalars describing the relative 
motion of the two D8 branes. Fortunately, we may easily add the to the 5D theory by 
applying the U{2) symmetry to the effective Lagrangian fl5.14f[5rT5|) . Thus, without re- 
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expanding the DBI action for two separated branes, we immediately obtain 



^Hn = ^^^^[{u{w)/uo)x[^^^^^^ (5.16) 
Ant = 4^ tr(($ = $ + ^V'^) X {-l{:Fp + 5{:Fp)) + 0($^) 

+ 2$^^ X + + o{^'). (5.i7) 

However, for the two-baryon system we are interested in, the isotriplet scalars fields are 
much weaker than the isosinglet field $ because they have much weaker sources. Indeed, 
for two static baryons, the SU{2) gauge fields are purely magnetic in 5D sense, i.e., Fq- = 
^ow = O5 while the U{1) gauge fields are purely electric, Fij = Fi^ = 0. Consequently, on the 
last line of eq. f l5.17p . both F^^^F^^^ = and F^^F^^ = 0, which leaves the isotriplet scalars 
with no source at all. 

For the baryons with non-zero spins, the SU (2) electric fields do not exactly vanish, but 
they are much weaker than their magnetic counterparts. Specifically, 

pSU{2) ^ p X 1/2 J .J 

( ^^p'° ] - ^ X — < 1 rs 18) 

FI^P V ^static/ Mbaryoni?baryon VX ' ^ ' ' 

At the same time, the abelian electric fields are also much weaker then the non-abelian 
magnetic fields. Indeed, were it not for the Chern-Simons interactions 

r - tr (.JKLMN A TT T -L \ 

= Q^^2 ^ ^''^^^^^ (^^{^klFmn) + \FklFmn^ 

= |A^c^j X (instanton current)"^ + (5.19) 

the baryons would not generate any abelian fields at all. As it is, for baryons of radius p, 
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the non-abelian magnetic fields are 



orr^nN I O [1 1 o^) ill tfic Heal zone, 

F^T = { ! (5-20) 

0{p/r ) in the intermediate zone, 

while the abelian electric fields are 

0{Nc/ KT"^) in the intermediate zone, 
where k = 0{NcXMj>^) is the 5D kinetic-energy coefficient, thus in both zones we have 

- - ' ' « 1. (5.22) 



As to the abelian magnetic fields, they are generated by the Chern-Simons terms involving 
Fmag'^^ X SO they are even weaker than the electric abelian fields. Altogether, we have 

a hierarchy of gauge fields 

FIT « ^"i''' « ^cf « F^il^ . (5.23) 

Consequently, the scalar-vector interaction Lagrangian fl5.17p provides a much stronger 
source for the isosinglet scalar field $ than to the isotriplet fields which leads to the 
scalar field hierarchy 

1 J , , 

— ~ - X — < 1. 5.24 
^ X Nc 

Hence, the nuclear forces due to the triplet are much smaller then the forces due to the 
singlet $, so we shall disregard the through the rest of this article. 

Focusing on the singlet scalar $, we see that the dominant source for it comes from the 
SU{2) magnetic fields, so to the leading order in 1/ A we may approximate 

Ant ^ X {-l^F^j)' + \{F^J') + unimportant. (5.25) 

Moreover, in the near and intermediate zones where the 5D gauge coupling is approximately 
constant — k{w) = kq x {u{w)/uq) ^ kq for \w\ <^ M^^, cf. eq. (15. 9p . — the SU{2) fields 
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are self-dual in the 4 space dimensions, 



Ftj{x,w) = eijkFU^,w) (5.26) 

Their specific form is given by the ADHM self-dual solution with instanton number = 2, but 
fortunately we don't need the gory details of this solution here. All by itself, the self-duality 
assures us that 



1{F;^Y = i(F;j'= 3^e°™^(F^^F^,^) ^ 47r^x instanton density /(x,^) (5.27) 



and hence, the scalar field $ couples to the same instanton number density as the abelian 
vector potential Aq, 

Ant = 1^ X (-3 X An^I + 5 X An^l) = ^ $ x /, (5.28) 

N - 

cf. £sc = -i^A^xI. (5.29) 

Therefore, in the near and intermediate zones where both the and the $ have approxi- 
mately constant kinetic terms 

CUKM - ^^^^ [\Fl^ + Y^T^ X 1(9.,$)^} , (5.30) 

both ^4° and $ are 5D Coulomb fields of the same charge density I{x,w), 

MS,w) = T^^Jd'^dw' ^, (5.31) 

NcXM\( J [x' — x)^ + [w' — wy 

$(x,ti;) = ^ X Aq{x,w). (5.32) 

This fixed ^/ Aq ratio has profound consequences for the attractive nuclear force: For any 
geometry of the instanton density I{x, w) in the near and intermediate zones, there is a fixed 
ratio between the attractive force due to $ and the repulsive force due to Aq, namely 

"^attractive C""^ 

-— r-. = —Ca/r = 7^ = COUSt. (5.33) 

^repulsive/' J, ^ "/^ ^2 ^ ^ 

Note that this ratio vanishes ioi ( = 1 — there is no attractive scalar force in the antipodal 
model. For the non-antipodal models, the attractive/repulsive force ratio Ca/r increases 
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with (, but it never gets larger than 1/9. Thus, in the near and intermediate zones of the 
SakaiSugimoto model, the attractive nuclear force is always weaker than the repulsive force. 



5.2 Attractive Forces in the Intermediate Zone 

In hght of eq. f l5.33p , calculating the attractive force between two nucleons at an intermediate 
distance from each other 

^-3/4_)^-l/2 ^-1/2 
P = ^baryon TT < r < W^ax T^J— (5.34) 

seems like a simple exercise. Approximating both baryons as point sources of 5D Coulomb 
fields Aq and we get 

1 — c^^ 37r/v n — c^-^) 

^attractive = ^ ^ yrepulsi.e ^ - ' \ , ^ ^ ' xl/r\ (5.36) 

and the net potential is a repulsive +l/r^ hard core. Note that both the repulsive and 
the attractive potentials (I5.35f[5l36|) are blind to spins and isospins of the two baryons. 
However, we saw in section 4 that such blindness is an artefact of the treating the two 
baryons as point sources. A better approximation makes the repulsive potential sensitive to 
the baryon's spins and isospins, and we shall see momentarily that the attractive potential 
has a similar sensitivity. 

Indeed, let's follow Hashimoto et al [16] and take the SU{2) gauge fields to be exactly 
self-dual, i.e., the two-instanton ADHM solution. The instanton density of this solution 

OKLMN , . 

l{xM = —-^tT(FLKFMN) = I^'\x,w) + I'^'\x,w) + AJ--"^P(f , w;) 

dln^ V / ADHM 

(5.37) 

differs from two separate instantons by AJ°™''''^p = 0{{p/rY), which has two effects: First, 
the repulsive Coulomb self-interaction of the instanton density (15.371) 

(5.3''^ 

is more complicated then fl5.35p . and second, for a w-dependent gauge coupling k{w 
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/t(0) X {u{w)/uq) — cf. eqs. fl5.16p and (15.91) — re-distribution of the SU{2) gauge fields in 
the w direction changes their energy by 

AE[SU{2)] = ^f^^J"^^ yd'xdw{u{w)/uo) X |tr(F,%) - 2 x 87r^| 

= ^!£^x^^^x [d'xdww'xAr^'''-^{x,w). (5.39) 
Sivr 3C J 



Both effects were evaluated in careful detail by Hashimoto et al [T6], so let us simply adapt 
their results to the present situation. For two nucleons — i.e., baryons of spin = isospin = ^ 
— and 1 (which suppresses the quantum fiuctuations of the baryons' sizes or locations), 

they found 

AE[SU{2)] = (^1 + y(li-l2) X [2(n- Ji)(n- J2) - (Ji-J2)]) (5.40) 

where A = x ^ x p^, (5.41) 

1627r 3C 

V[U{1)] = (^Uaivc - ^ + §(li-l2) X (Ji-Js)^ (5.42) 

Or rather, they obtained these formulae for the antipodal modal, but the extra C^dependent 
factors for the non-antipodal models are obvious from eqs. f l5.39p and (15.381) . 

Now consider the attractive scalar force between two baryons. Since the scalar field 
$ couples to the same non-trivial instanton density (15.371) as the abelian gauge field Aq, 
the attractive force has the same complicated spin and isospin dependence as the U{1) 
force (I5.42p : Instead of the naive eq. (15.360 , we get 

Vm = -C,/rXV[U{l)] = -^^x (^Ui^e - I + §(li-l2) X (Ji-J^)) (5.44) 

where Ca/r = |(1 ~ C~^) as in eq. (I5.33p . 

But besides the direct contribution fl5.44p of the scalar force to the two-baryon potential, 
it also affects the baryon radius p — which in turn affects the SU{2) force according to 
eqs. (15.40115^^1]) . To see how this works, consider a stand-alone SU{2) instanton with 

I{x,w) = — ^- — (5.45) 
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The classical energy of this instanton is 



E{p) = E[SU{2)] + AE[U{1)] + AEm (5.46) 

where 



(5.47) 

the Ao-mediated Coulomb self-interaction is 



AE[f/(l)] = +fi X (5.4J 
5p"^ 



where i? is as in eq. f l5.43p . and the ^-mediated Coulomb self-interaction is 

AE[$] = -C,/, X AE[[/(1)] = -C,/,5x A. (5.49) 

Note that the scalar / abelian vector force ratio f l5.33p works at short distances as well as 
intermediate. 

Assembling all the contributions to the classical instanton energy, we get 

and minimizing this formula with respect to p gives us the classical mass and radius of the 
baryon: 

M,.,,o. - ^7^^ + 71 l,^C^ J x(l-a/.) , (5.52) 

Note that these formulae differ from eqs. (13.181) by factors (1 — Ca/rY"™'^'^"'"'^^'- These factors 
are due to ^-mediated attractive force which we didn't take into account back in section 3. 
In particular, the scalar force reduces the baryon radius by a factor {\ — Ca/rY^^- Substituting 
this radius into eq. (I5.4ip gives us 

v4 = [l-Ca/r), (5.53) 



i?baryon = P@mm = | ^J^^ 1 xj^^^^— ^) X (l - Ca/,) , (5.51) 
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which means that the indirect effect of the scalar force reduces the S'f/(2)-mediated nuclear 
force by the factor (1 — Cq/^.). At the same time, the direct effect fl5.44p of the scalar force 
reduces the net Coulomb force by exactly the same factor, 

V[U{1)] + = V[U{l)]^{l-Ca/r). (5.54) 

Altogether, we end up with this (1 — Ca/r) factor multiplying the whole nuclear force in all 
its complicated glory, 

V;et(r) = AE[Sf/(2)] + V[U{1)] + 



^x(l-a/.)x (^1 + ^(Ii.l2)x [2(n.jO(n-J2) - (Ji-J^)]) 
B / 2 32 \ 

+ -2^{^' C-h) X (Inaive - - + - (ll ■ I2) X ( ■ J^) j 
^X (1 - Ca/r) X (^Inaive + ^ (ll " I2) X (n • Ji) (n ■ J2) ^ (5.55) 



To summarize, the net holographic nuclear force in the intermediate zone is precisely as 
in Hashimoto, Sakai, and Sugimoto [16]. The only effect of the attractive force due to the 
scalar field $ is to reduce the whole force by a constant overall factor (1 — Ca/r)- Depending 
on the C parameter of the model, this factor varies between | (reduction by 11%) and 1 (no 
reduction at all), but it is always positive so the net force in the intermediate zone is always 
repulsive. 



5.3 Attractive Forces in the Far Zone 

In the far zone, the if-dependence of the gauge coupling 

^ , ,5,6) 

2167r'^ Mo 

becomes important and the lowest-energy two-instanton solution of the SU{2) equation 
Dm{K'F'^^) = is no longer self-dual. Consequently, working out the overlap between 
two distant instantons becomes rather difficult and we are reduced to a cruder point-source 
approximation. Earlier, we saw that in the intermediate zone, this approximation has yielded 
a correct isoscalar central force but missed the isovector spin-spin and tensor forces. By 
analogy, we expect that in the far zone, the point-source approximation would give us a 
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correct isoscalar force, but the isovector forces could be wrong. 

Let's focus on the isoscalar forces. In the point-source approximation, they follow the 
Green's functions of the Aq{x,w) and w) fields in the far zone. 



yrepulsivc^^^ 



+ 



X 



CAMa 
X (x, 



(5.57) 



uiw 



Uq 



V' + 



d ( u{yS)\ d 



dw \ Uq J dw 



1 -1 



0,0,, 



V 



attractive 



(r) = -Ca/r X 



X 



(5.5^ 



X (f, 



u(w . 



Uq 



d f u{w)\^ d 



dw \ Uq J dw 



+ ^i\w) 



0,0 



where 



2CM1 
9 



X 
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u[w , 



Uq 



10 



32C" 



u[w] 



Uq 



+ 13(1 -C 



U[W] 



Uq 



(5.59) 



is the effective mass term for the $ in the expansion of the DBI Lagrangian in terms of 5D 
coordinate {x^, w)E Expanding the 5D fields Aq{x, w) and $(x, w) in terms of the 4D vector 
and scalar mesons, we may express their Green's functions — and hence the potentials f l5.57f - 
I5.58P — as sums of Yukawa potentials 



^repulsive -J 
■^attractive 



' 2CAMa 



Y^\K{w = ^)\'^'-^^t<ll 



n=l 



El ^ c , , |2 expl— m„r) 

|^n(^ = 0)| X " ^ 



(5.61) 
(5.62) 



n=l 



Here \E'Jf (w) and '^^{w) are the wave functions in the fifth dimension of the respective vector 
or scalar mesons. Note that only the odd-ra modes — which give rise to the true vector 
1~ and true scalar O"*" mesons — contribute to the sums (l5.61f[5T62l) . The even-ra modes — 



®The negative sign of /i^(w) is an artefact of the w-dependent normalization of the scalar field; a canon- 
ically normalized $ has positive mass^ 



w) 



CM 



20 



u(w) 

Uq 



Uq 



44(1 - C"^) 



u{w) 

Uq 



(5.60) 
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responsible for the axial-vector l"*" and pseudo-scalar mesons — don't contribute because 
their wave functions vanish aX w = 0. 

We do not have analytical formulae for the meson's masses or wave functions, but numer- 
ical calculations [HI |32] show that for Sakai-Sugimoto models with any C ^ I5 the vector 
mesons are always lighter than the scalar mesons, 

Vn, m]^ < mf. (5.63) 

Numerical calculation of the mesons' wave functions at w = is still in progress, but there 
does not seem to be much difference between vector and scalar mesons, so it is reasonable 
to assume that 

|^n(0)|^ 

for odd n, - — ^ > Ca/r = • (5.64) 

1^^(0)1' "^'9 ^ ^ 



If this assumption is correct, then every term in the repulsive potential f l5.6ip is stronger 
than the corresponding term in the attractive potential (I5.62p . and the net isoscalar force is 
repulsive throughout the far zone. 

At the outer end of the far zone, we don't need the assumption (I5.64p to show that the 
repulsion is stronger than attraction, all we need to know is that the lightest vector meson is 
lighter than the lightest scalar meson, mX < mf. Indeed, for r ^ (l/'n^meson) ~ (1/Ma.a/C), 
each sum (I5.61f[5l62l) is dominated by the slowest-decaying Yukawa term belonging to the 
lightest vector or scalar meson, 

^ ^ (5.65) 

Regardless of the pre-exponential factors here, the longer-ranged force always wins over the 
shorter-ranged force at long distances, and since < mf in all versions of the Sakai- 
Sugimoto model — antipodal and non-antipodal with any C > 1 — the net isoscalar force 
is repulsive at long distances r ^ 1/ (Ma-\/C)- 

Note that this behavior of the Sakai-Sugimoto model is very different from the real-life 
nuclear physics. Indeed, in reality the lightest scalar meson a is lighter then the lightest 
vector meson, ~ 600 MeV while ^ 787 MeV, and consequently the real-life isoscalar 
nuclear force is attractive rather than repulsive at short distances. We do not know the reason 
for this discrepancy; perhaps it's a peculiar bad feature of the Sakai-Sugimoto scheme and 
could be fixed by an alternative holographic model. But perhaps it's an inherent problem 
of the large-iVc nuclear physics. As explained in section 2, the QCD origin of the a (600) 
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resonance is rather controversial — maybe it's a true meson originating in the a field of the 
linear-sigma-model-like chiral symmetry breaking, or maybe it's just a two-pion resonance 
which would not exist without strong vrvr interactions. In the first scenario, the large Nc limit 
would make cr(600) into a narrow resonance, but it would still be there and lighter than the 
lightest vector meson so the attractive nuclear force would have a longer range than the 
repulsive force. But in the second scenario, there would be no cr(600) resonance in the large 
Nc limit, the lightest remaining scalar mesoij^ would be heavier than w, and the dominant 
isoscalar nuclear force at long distances would be repulsive rather than attractive. 

The best way to settle this issue would be to find the a resonance and its mass in a lattice 
QCD calculation for several values of N^-, which would hopefully allow us to extrapolate to 
Nc ^ oo. Alternatively, once we have several different holographic models, we can compare 
their predictions for the meson spectra in general and for the lightest scalar to lightest vector 
mass ratio in particular. Either way, this issue will have to wait for future research. 

Another issue that would have to wait for future research involves nuclear forces arising 
from nucleons exchanging pairs of un-bound mesons — especially pairs of pions — rather 
than single mesons or resonances. In the real-life nuclear physics, the double pion exchange 
generates the longest-ranging attractive isoscalar force, which makes a significant contribu- 
tion to the bulk binding energy of the nuclei. In the large Nc limit, the double-pion exchange 
decreases as 1/Nc relatively to the single- meson exchanges, but it would remain significant 
at longer distances where the exp(— mr) factors for single mesons like a or u are even smaller 
than l/Nc- Unfortunately, in holographic duals of QCD, the double-meson exchanges hap- 
pen at the one-string-loop level, so they cannot be be calculated in terms of an effective 
semi-classical 5D gravity — or even lOD gravity plus other local fields — but require a fully 
quantum string theory. Since string perturbation theory in curved backgrounds is rather 
hard, we leave the double-meson exchanges for future research. 

Instead, lets us now address the dominant nuclear force at the longest distances, namely 
the isovector force due to single pion exchange between two nucleons. 



J1J2 X -i^ + Ti2(n) X + ^ + - 

Sr \ 6r J 



> -Jl^ (Iil2)ri2(n) X - (5.66) 



^ T 



^Probably the /o(980), or maybe even the /o(1450) — which is the best fit for the lightest scalar in the 
Sakai-Sugimoto model — in case the /o(980) is also a two-pion resonance that would go away when Nc ^ co. 
Either way, the lightest surviving scalar would be heavier than w(787). 
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where Ti2(n) = 3(nJi)(nJ2) — (J1J2) is the direction dependence of the tensor force. The 
overall coefficient of this force was calculated by Hashimoto et al 116] for the antipodal 
model. Adapting their method to non-antipodal models gives us 



7r/2 36vrF(C) 



' X (5.67) 



where p is the classical radius of the baryon and F is the normalization factor of the 
pion's wave function 

<i/^{w) = , ^ X (5.68) 

m - ^ f'l^-^- Uvi ^ ' , (5.69) 



3 1 u(v!) J "Vy'-C-V-li-C 

~W?niax 1 

The pions themselves are modes of the SU{2) gauge fields A1.j{x,w) and don't know the 5D 
scalar field $(x,w) from the hole in the ground. However, their interactions with baryons 
depend on the classical baryon radius p, and because of the ^-mediated attractive forces in 
the near zone, this radius is smaller than it would have been otherwise, 



8l7rv/2/15\ 



2 



Consequently, the pion-mediated nuclear force in the far zone becomes sensitive to the $- 
mediated attractive force in the near zone, 

_ (Iil2)Ti2(n) 48iV. <3/2 
^ ^"^^ - r3 5AM1 (8C3-5)F(C)^ ^"Z''^- ^^'^^ 

What about the other isovector mesons' contributions to the nuclear force? In the 
Sakai-Sugimoto model we have (pseudo) scalar isovector mesons coming from the modes 
of the 5D scalar fields ^"'{x,w), but their couplings to baryons are too small to create 
an appreciable force. In addition, we have vector and axial vector isovector mesons p 
(770 MeV), ai (1260 MeV), etc., coming from the non-zero modes of the SU{2) gauge fields 
A1j{x,w). Their contributions to the isovector spin-spin and tensor forces were calculated 
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by Hashimoto et al [16] as 



1 Ofi 3 °° 



27 

71 = 1 



(J.J.) ^ + r,.(„) (4 + + 

(5.72) 



where 



+ |^„(w = 0)r for odd n, 
|^n(w^ = 0)1 for even n. 

Note the overall factor for all the isovector forces. Since the baryon radius is affected by 
the attractive forces in the near zone, they have indirect effect on all the isovector forces in 
the far zone, 

^isovector OC (l-Qr). (5.74) 

This rule applies even to the effects of the baryon-baryon overlap that Hashimoto et al could 
not calculate in the far zone. We cannot calculate their r dependence either, but their p 
dependence should be the same as in the intermediate zone, 

V^overlap^^^ OC OC [l-Ca/r)- (5.75) 



6 Full DBI Action in the Near Zone 

In the previous section we studied nuclear forces in the intermediate and far zones, but since 
those forces depend on the baryon radius p, we had to stick our noses into the near zone to 
see how p is affected by the scalar-mediated forces. In the intermediate and far zones, all 
the gauge fields are weak enough so we could expand the DBI Lagrangian in power of J^mn 
and truncate the expansion after the leading Yang-Mills term as we did in eqs. (]5.13f[5n5l) . 

for r > p ~ —= , ItkJTmn ^ Gmn 

VAMa 

I , , ^'-'^ 

£dbi oc str Vdet(^ + lixa'T) ^ ^- det(^) x + ^^^^ ^ \,y [TmnJ=^^^)^ ■ 
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However, in the near zone r p the SU{2) gauge fields become too strong for the YM 
approximation, 



for r ~ p 2na'TMN ~ 9m n =^ Cbbi ^ v^^dit(^ (^2 + ^^^^tr {TmnJ^'^'')^ , 

(6.2) 

which casts doubt on accuracy of eqs. f l5.5H - [532l) for the baryon's mass and radius and hence 
of all the p-dependent formulae for isovector nuclear forces. 

In this section we shall see that despite the higher-order DBI interactions beyond the 
Yang-Mills approximation, to leading order in 1/A the baryon radius remains exactly as in 
eq. (I5.5ip and we do not need to make any leading-order corrections to the nuclear forces we 
have computed in the previous section 4. And at this happy note, the readers impatient with 
technical details may skip the rest of this section and go straight to the summary section 7. 

Our main point is that while the Yang-Mills approximation (16. ip to the full DBI actioi|l° 
does not work for generic strong gauge fields, it may work 'by accident' for some special 
J^mn{x) configurations. In particular, for self-dual gauge fields living on a fiat D4 brane 
stack, the YM approximation not only works but happens to be exact, regardless of the 
number of instantons or their sizes |35j , 

for J-^, = |e^,„/3^"^ str Vdet(l4x4 + 27ra'^) = 2 + tr ( J^m^^*^^) (6.3) 

In our case, the 5D metric is not fiat, and besides the (approximately) self-dual SU{2) 
magnetic fields we also have the abelian electric field. Nevertheless, eq. (16. 3 p continues to 
hold exactly for a slightly modified self-duality condition for the SU{2) magnetic fields. 

Indeed, consider a more general case of some non-fiat but static metric 

ds^ = -\goo{x)\dt^ + gmn{x)dx"^dx'^ (m, n = 1, 2, 3, 4), (6.4) 

arbitrary but purely-electric abelian fields Fmo{x) and purely-magnetic SU{2) fields F^^{x). 
In general, Tseytlin's non-abelian version [34J of the DBI Lagrangian works like this: First, 
one calculates the determinant of the gMN + 21^ a' Tun matrix (in space-time indices) while 
completely ignoring their gauge indices or the fact that they don't commute with each other. 

""^^ Actually, the DBI + CS action is also incomplete — for multiple D-branes, there are additional terms 
involving covariant derivatives 'Dj^Tmn of the non-abelian gauge fields. Fortunately, for instanton-like gauge 
fields of a baryon, such covariant derivatives are relatively small — even when the gauge fields themselves are 
large as in eq. (|6.2p — so following Tseytlin [33] we shall limit our analysis to the non-abelian DBI action. 
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Second, one expands the square root of this determinant into a formal power series in the 
gauge fields. Finally, for each term in this determinant one takes a symmetrized trace over 
the gauge indices, and then tries to re-sum the series. For the static case at hand we have 



MN 



I - 


900 




-Tca'Fmo \ 




9mn ~l~ 27rQ; Frnn J 



and hence 



where 



det(5fA/iv + 27ra'J'A/iv) = Is'ool x det(^rnn + 27ra'Fi^n) 

5D ^ ' 4D ^ ' 



9mn{x) = 9mn{x) - (vTa') \g | F^oi^nO • 



(6.5) 

(6.6) 
(6.7) 



Note that while the determinants on both sides of eq. (16. 6 p are formal — they ignore the 
non-commutativity of the magnetic fields Fmn = \'t°'F^^ — but the modified metric (16. 7p 
does not have any non-commutativity problems because the electric fields F^q are purely 
abelian. When those electric fields are too strong, we can get a different problem of Qmn 
matrix loosing positive-definiteness, but fortunately the U{1) fields of a baryon never get 
that strong: Even in the near zone. 



{Tia'Y\d''\g^''F^,Fn, < 0(1/A) « 1, 

so the modified metric Qmn remains safely positive-definite. 

On the right hand side of eq. (16. 6p . the 4D determinant evaluates to 



(6.8) 



det{gmn + '27!-a'Fmn) = det{gr, 
where 



1 + 2i7ra'frTF„,nFp, + (vra')' [r^r^FmuF, 



pq 



Fpq = l^/d^ep.rsf'^r-Frr 



(6.9) 

(6.10) 



is the Hodge dual of the Fmn with respect to the modihed metric (jmn- When the SU{2) 
magnetic field is self-dual (or anti-self-dual) with respect to that metric, F^n = ±-^mn, the 
4D determinant (16. 9 p becomes a full square, hence 



det{gMN + 2na'J^MN) = v^Vdet(^) x 1 + (vraO^^'^P^'i^^nF; 

5D L 



(6.11) 
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and the symmetrized trace becomes the ordinary matrix trace (over the U{2) gauge indices). 
Consequently, the complete DBI action for all the gauge and metric fields splits into the 
ordinary Yang-Mills action for the (self-dual or anti-self-dual) SU{2) magnetic fields, plus 
the abelian DBI action for the metric and the U{1) electric fields only, 

Cdbi = -T8e-<^Vol(^^)v^x^^(±6™)F^„F« 

-2T8e-'^Vol(5^) X y-det((7Mjv + vra'FMiv) (6.12) 

(where on the second line we have used |(7oo| x det4D (fi'mn) = — ^Q^bT:>{.gMN + t^O-'Fmn))- 

For the baryon, the SU{2) magnetic fields become strong in the near zone, but the U{1) 
electric fields they induce (via the Chern-Simons interactions) are relatively weak {0(1 /\/X) 
or weaker) in all zones, and the perturbations of the 5D metric due to the scalar field $(x) 
are also weak in all zones. Consequently, the second line of the Lagrangian fl6.12p gives rise 
to the usual kinetic terms for the scalar and the abelian gauge fields, while all the higher- 
order terms are suppressed by negative powers of the 't Hooft coupling A (and outside the 
near zone also by positive powers of (p/r) ^ 1). At the same time, the first line of the 
Lagrangian fl6.12p gives rise to the Yang-Mills Lagrangian for the SU{2) gauge fields, and 
also to their interactions with the scalar field $ and with the background metrics's curvature 
via the w-dependence of the 5D gauge coupling 

k{x,w) = 2{narTse~t'Vo\{S')^\ = + ^x,w)] . (6.13) 

However, we do not get any terms with higher powers of the strong SU{2) gauge fields, and 
everything works precisely as in §5 — except that the F^„(x, w) fields should be self dual 
with respect to the modified metric Qmn rather than the true metric Qmn- 

To see the effect of this modification of the self-duality condition, consider an SO {4) 
spherically symmetric instanton-like field configuration 

Aq{x,w) = 0, Am{x,w) = h{r) x —iU'^dmU, U{x,w) = , r = \/xP + w'^ 

r 

(6.14) 

where m = 1, 2, 3, w and h{r) is some smooth function of the 4D radius. The corresponding 
magnetic field strength has form 

Fmn{x,w) = -id^^h X U^dn]U + lh{l - h) X [U^ d^U, dnU] . (6.15) 
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In particular, along the w axis we have 



F". = -F"" = -—xS" F". = ^ ^ X e"'^ (6 16) 

^ wt ^ tw r dr ^ ' ^u.iu; 



while in other directions we have similar fields rotated by appropriate SO {4) x SU{2) sym- 
metries. The instanton density of such fields is 

J ^ 6^^-F,",F^„ ^ 3 hjl -h) dh ^ ^g^^^^ 
647r^ TT^ dr ' 

so the net instanton number is 

oo fe(oo)=l 

#Instantons = J I (r) x 27r^ dr = j 6h{l - h) dh = 1 (6.18) 

h{0)=0 

for any radial profile h{r) that satisfies the boundary conditions. 

The magnetic fields fl6.15p become self-dual = \emnpqF!^q with respect to the Hat 4D 
metric Smn when 

r dr 

solving this differential equation gives us the usual instanton profile 

/iinstanton('") = ^- — n ioT some constant p. (6.20) 

_|_ pZ 

The self-duality condition with respect to the modified 4D metric calls for a slightly 
different profile 

+ X (1 + c)(r)) 

where the correction 6{r) depends on the abelian and the scalar fields. We shall see mo- 
mentarily that this correction is rather small, 6{r) = 0(1/A), so the baryon profile (16.211) is 
approximately the usual instanton profile (I6.20p . 

To obtain the modified metric Qmn and the corresponding self-duality condition we need 
the ordinary 5D metric Qmn- Eq- (IS.lOp gives us this metric to first order in the scalar field 
$, and that's a good enough approximation even in the near zone where 

<!>{x,w) = 0(Ma) < ^ = O(AMa). (6.22) 

a' 
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Consequently, expanding the modified metric 



in powers of 1/A (in tlie near zone wliere 6u{w) = u{w) — Uq < (uq/X)), we get 

3/2 



(6.23) 



ds' 



R 



Di 



Uo 



R 



DA 



3/2 



dx"^ + dw'^ + 0(1/A) 
36u + 37r«'$ 



1 



+ [Tca' 



2uo 

Uo 



,^2 A 35u~ 13na'<i>\ , ^ 

dx^ + (1 A dw^ 

V 2mo J 

(^YZ^l^-*^^^")' - {dn^Aodx^Y^ 



+ 0(1/A2) 



(6.24a) 
(6.24b) 
(6.24c) 
(6.24d) 



Note tliat to tlie zerotli order in 1/A, tlie modified metric (16. 24b ) of tlie near zone is simply 
flat, so the modified self-duality condition is just the good old flat-space self-duality condition 
^mn = \^mnpqFpq. And that's why the DBI correction 5{r) to the baryon profile (I6.2ip is 
0(1/A), hence 



i?baryon[full DBI] = /?baryon[aS lu Cq. ESI] X (1 + 0(1/A)), 



(6.25) 



and all our results of §5 for the intermediate-zone and far-zone nuclear forces are indeed 
correct to leading order in 1/A. 

For completeness sake, let's calculate the modified baryon profile (I6.2ip to first order in 
1/A using the modified metric fl6.24b -d). Since the overall conformal factor of gnm does not 
affect the self-duality condition (16.101) . let's factor it out as 



C 



uiw , 



R 



D4 



3/2 



27r«'$ 



u 



(6.26) 



and focus on the remaining deviations of the cjmn from flatness 

2W$ 



ds' 



where 



C X 



{dx^ + dw^) + V{r) X dr^ + 



u 



X [dx^ -3dw^) + 0(1/A2) 



Vir) 



TTO; 



, i\2 



Mo 



1-C 



-3 



- {dAo/dr) 



(6.27) 



(6.2^ 
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Note that the first two terms inside the square brackets in fl6.27p are spherically symmetric 
in all 4 space dimensions, but this 5*0(4) symmetry between the x and w coordinates is 
broken by the third term. Consequently, a single static baryon is SO {4) symmetric only 
to the leading order of the 1/A expansion but the sub-leading terms spoil this symmetry. 
To see this asymmetry we should use a more complicated ansatz than (16.141) for the non- 
abelian gauge fields, so let's leave this issue for future research. For now let us focus on the 
baryon's radial profile, and in first-order perturbation theory this profile depends only on 
the spherically symmetric part 



symm 



C X (^(1 + V{r)) X dr^ + x dQ^^ (6.29) 



of the modified metric gmn- 

For a spherically symmetric metric ( 16.29P the self-duality condition f l6.10p becomes very 
simple: The SU{2) magnetic fields with one radial and one tangential index should be 
a/I + V times stronger than the fields with two tangential indices, for example along the w 
axis we should have 

= VlTV X le,j,F^, . (6.30) 

In terms of the radial profile h{r) of the spherically symmetric SU(2) fields fl6.14p . the 
modified self-duality condition amounts to 

2 dh I — - 4/i(l -h) , , 

__ = /rTPMx^^. (6,31) 

Solutions of this differential equation have general form 

h{r) 



/2dr 
— ^Jl + V{r) + const, (6.32) 



1 - h{r) 

and for Vij-) = 0(1/A) ^ 1 as in eq. (I6.28p we may approximate these solutions as 



^2 



+ p^(l + o(r)) 

where 

oo 

rV(r')dr' 

6{r) = j ^ j . (6.33) 

r 

In the first-order perturbation theory we may calculate D(r) using the scalar and electric 



45 



fields of tlie un-perturbed baryon. In tlie near zone 



lience 



$(r) = ^ X Ao(r), (6.35) 



J \XMaJ (r2 + p2)6 



TT 



(40C=' - 25)3/2 (r/p)2(3 + (r/pff 



and 



A 23/2C5/r^I^ 4(1 + (r/p)2)6 (6-36) 

_ - (40C^ - 25)3/2 32 + 25(r/p)^ + 5(r/p)^ 
^ ^ " A 23/2CV1-CV 120(1 + (r/p)2)5 • ^6-37^ 

As promised, this correction to the baryon profile is small in the near zone — 5(r ~ p) = 
0(1/A) — and becomes even smaller at larger radii. 



Thus far, we have merely assumed the self-duality of the SU{2) magnetic fields with 
respect to the modified 4D metric gmn- To justify this assumption, we are now going to show 
that the lowest-energy field configuration with A^instanton = 1 is indeed self-dual, or at least 
approximately self-dual to leading order in 1/A. 

For simplicity, we look for the minimum of energy among the 5*0(4) symmetric field 
configurations only, i. e. we restrict the SU{2) gauge fields to the ansatz f l6.14p but allow for 
generic radial profiles h{r) (subject to boundary conditions h{0) = 0, h{oo) = 1). Likewise, 
we take the modified 4D metric Qmn to be spherically symmetric as in eq. (16. 29 p . For such 
configurations, the 4D determinant (16. 9 p becomes 



det [g 

mn + 27Ta'Fmn) = det{g 

n 



^ 2 

pq 



= C^(l + I?) X (l + aV')(l + /3V2) (6.38) 

where 

27Ta' 2h(l — h) ^ 27ra' 1 dh , „x i /o 

« = -^x^-^, /3 = -^x--x(l + P)-V2, (6.39) 

and = rf + T2 + t"^. Please note that we are not allowed to use the Pauli matrix algebra 
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and set = 3 while we calculate the determinant fl6.38p : instead, we should to treat the Tj 
as if they were independent and un-constrained commuting numbers. 

Given the determinant fl6.38p . we need to expand its square root in powers of arj and 
/Sxj, then for each term in this expansion we should restore the SU{2) gauge indices of the 
Ti,2,3 factors and calculate the symmetrized trace, and then we need to re-sum the series. 
This is easy to do in the self-dual case of a = P when the determinant (16.381) is a full square, 
but for general a ^ (3 this calculation takes a few pages. We present it in the Appendix to 
this article; here is the end result: 



str V(l + a2f2)(l + /32f2 = . (6.40) 

For a given product ax/? this expression is minimized for self-dual a = P — this is the 
DBI version of the BPS lower bound on the Young-Mills action for gauge fields of a given 
instanton number. For convenience, we may rewrite eq. (I6.40p as 



str ^/{l + a2f2)(l + /32f2) = 2 + 6af3 + P{a, f3) x {a - Pf (6.41) 

where P(a, /3) is some complicated expression. It's gory details will not be important in the 
following, expect for the special case of approximately self-dual a ^ (i when 

Pia ^ /3) ^ + . (6.42) 

^ ^' (l + a2)(i + /32) ^" ) 

More generally, P(a, /3) is always positive and never greater then 3. For weak fields a, /3 ^ 1, 
P ^ 3 — indeed, the Yang-Mills Lagrangian tr(l + \{a + Pfr^) = 2 + 3(a + /3)^ has 
form (I6.4ip for P = 3 — so the difference 3 — P measures the importance of the higher-order 
DBI corrections beyond the Yang-Mills approximation. For example, the approximately 
self-dual baryon with h ^ r^/ (r^ + p^) fQ^- p eq. (|5.5ip has 

27ra' 2p2 \A(C) x . w.x 40(^-25 , 

a(r) « 0(r) « — X = for MO = (6.43) 

and hence 

5k 2A;2 
At the very center of the baryon, P(r = 0) dips to 2.20 for the antipodal model — and 
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even lower to 1.37 for the non-antipodal models with C ^ 1 — which indicates fairly strong 
higher-order DBI corrections. However, just one unit of p outside the center, P(r = p) 
climbs back to 2.8 for the ^ 1 models and even higher for the antipodal model, so the 
higher-order DBI corrections are important only in the inner core r < p of the baryon's near 
zone. Outside this inner core — from the outer side of the near zone all the way to the 
far zone — the higher-order DBI terms are small and one may safely use the Yang-Mills 
approximation to the SU{2) fields' Lagrangian. 

Plugging eq. (I6.4ip into the non-abelian DBI Lagrangian, we obtain 



-^DBi = T8e-'^Vol(5^) X str ^det{gmn + 27ra'F„ 

TgVoK^^ 



TgVoK^^ 



CVl + 1? X (^2 + 6a/3 + P(a, (3){a - (3f 



(6.45a) 
(6.45b) 



X 



, 2 _ ^(1 — ^) dh 



2CVl + 1? 
+ 12(27ra')' x 

+ {2.ar^^ X f i ^ 



2h{l - h) 



Vl + V 



(6.45c) 



Inside the brackets on the bottom line fl6.45cp of this formula, the first term gives rise to the 
DBI Lagrangian for the scalar and abelian vector fields but does not affect the SU (2) gauge 
fields. The other two terms generate the Hamiltonian for the static SU (2) magnetic fields in 
a fixed background of the other fields $(r) and Ao{r). The second term provides the minimal 
BPS energy of the self-dual fields, while the third term is the energy cost of deviations from 
self-duality: It vanishes precisely when the modified self-duality equation f l6.3ip is satisfied. 
Adding the Chern-Simons term 



n 



cs 



An(r) X I(r) 



3N^Ao{r) h{l - h) dh 

X ~ 



27r^ 



dr 



(6.46) 



to the Hamiltonian for magnetic SU(2) fields in a fixed background of Aq and $ fields, we 
obtain 



n 



/i-independent terms 



u 



2tt'^ \TTa' 



+ $ + 3 A 



X 



h{l 



247r2 Vvra' 



VT+V 



- h) dh 
dr 
1 dh 

r dr 



2h{l - h) 



X 



VT+v 



(6.47) 



48 



Using spherically averaged brane geometry 



u{x) = uq + 5u{x), 

6u{w ^ r) ^ ~ Xa'Ml x (w^ ^) (6.48) 

81C V 4 y 




^ X (1 - C./,) X (r/pf 



10(2 

in the near zone, we integrate this Hamiltonian density to 

H[h(r)] = /i-independent + — ^— ^ [ dr h(l — h)^ (6.49a) 

7ra' Jq dr 

oo 

+ nJ drh{l-h)'^x ( + 3io(r) + $(r)^ (6.49b) 
J dr \ na' J 



oo 



+ I i„Pia(rmr)) , (dh _ ^^^^ _ ^ yTTP^ ^ ,.„(! + 0(1/A)) 







(6.49c) 



Note that the integral on the top line (16. 49b ) has the same value | for any baryon pro- 
file h{r) satisfying the boundary conditions; this integral provides the leading contribution 
{\NcM\C) /27t\: to the baryon's mass. The second line (16.49b ) corrects the BPS energy of 
self-dual SU{2) fields in a non- uniform background of Ao(r) and $(r) fields, and also brane 
curvature 5u{r). This extra energy has a non-trivial dependence on the baryon's radius and 
profile; its overall mag nitude is 0{NcMi^) = O(A-^Mbaryon)- Finally, the third line fICTb ) 
is the energy cost of deviation from self-duality. The w-dependent factor here is Uq (rather 
than much smaller 6u + 7ra'(3Ao + $) on the second line), so the energy cost of a major 
non-self-duality would be O(Mbaryon)- And minimizing this extra energy is precisely why the 
baryon profile is approximately self-dual. 

Indeed, minimizing the Hamiltonian (I6.49P as a functional of the instanton profile h{r) 
gives us a rather messy differential equation 



d 2(1 -2/i) , 

— + ^ -Vi + v 

dr r 



r dh 

P{a{rlP{r))^[-j===--2h{l-h) 



(6.50) 



d (5m + 7ra'(3io + $) 

-Qh[l — ft) X ^ 

dr Uq 
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In the near zone, the ratio {Su + ■ ■ ■)/uq on the right hand side is 0(1/A) small, so on the 
left hand side we should have a similarly small 



dh 



y/TTV dr 



2h{l-h) = 0(1/A) 



(6.51) 



deviation from self-duality with respect to the modified metric c/mn- The net energy cost of 
this deviation is 0(A~^Mbaryon), which is too small to be concerned with at the present level 
of analysis. Likewise, the effects of self-duality violation on the baryon's average radius — 
and hence on the nuclear forces in the intermediate and far zones — are minor corrections of 
relative order 0(1/A) to the leading-order effects we have calculated in the previous section 5. 

To conclude this section, we notice that while the deviation of the baryon' profile from 
the modified self-duality condition is 0(1/ A) small, it is comparable to the modification of 
flat-space self-duality condition due to V{r) ^ 0. So for completeness sake, we would like 
to calculate both effects on the baryon profile to the same order 0(1/A). Let's parametrize 
deviations of /i(r) from a flat-space instanton using bir) as in 



h{r) 



r'^ + X (1 + 5(r)) 



Substituting this formula into the differential equation (I6.5UI) and expanding to first order 
in 5{r) and P(r), we get 



d_ 2{p^-r^] 
dr r{p^ + r^) 



(p2 _l_ nr'^Y y 



6p r 



(p2 + r'^y 



X 



dr 



Uq 



(6.52) 



where on the left hand side P{r) is as in eqs. f l6.43f[6l^^ and on the right hand side 



Su + 7r«'(3io + ^) 



T 



Uo 



2 + (p2 + ^2)2 



5P 



where T 



9n /AOC - 25 



- Ca/r = 0(1/A), 



(6.53) 
(6.54) 



Solving the differential equation fl6.52p is a straightforward exercise in calculus. Integrating 
the outer differential operators gives us 



r— + V 

dr 



T 2r\3r^ + 12rV + 14p' 



P r) 



X 



5p2(r2 + p^'i^ 



(6.55) 
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Note that the right hand side here is of the same magnitude as the T> term on the right hand 
side which modifies the self-duahty condition; indeed, in present notations 



V 



-2kT X 



(r2 + p2)6 



(6.56) 



Altogether, 

Sir) = 



T I J^L±^(kp\'ip' + r^^'' 



2 I 



[r^ + p 



5p2[3(p2 + r2)4 + A;p8] 



(6.57) 

where the analytic form of the integral is rather unwieldy. Instead of writing it as a formula, 
let's plot 6{r) for 3 representative values of (: 




(6.58) 



Blue line for C 



1.283. A' 



red line for ( ^ 1, A; = |; black dots for the YM limit Ii6. 59\) . 

The differences between the three colored lines here are due to the higher-order (beyond YM) 
terms in the DBI action whose effect depends on the ( parameter. The low-order interactions 
— Yang-Mills, Chern-Simons, ^J-'^, and the D8-brane curvature 6u{r) — have the same 
radial dependence (in units of p) for all the Sakai-Sugimoto models, so when the higher-order 
interactions pucker out at r > p, the deviations from self-duality become ^independent and 
all the colored lines converge to the same black-dotted line 



T / r2 
Mr) « ^ 



P- 



J.2 _|_ p2 



(6.59) 
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This line shows that the low-order interactions themselves make the SU(2) fields deviate 
form self duality. In particular, at large radii we have a growing deviation 6 oc — due to 
brane curvature 6u oc (and hence the 5D gauge coupling growing with r). Consequently, 
self-duality of the SU{2) fields becomes less accurate in the intermediate zone r ^ p, and 



7 Summary and Open Questions 

Viable holographic nuclear physics obviously requires both attractive and repulsive nuclear 
forces. In holographic context, the hard core repulsive potential was found in [16], and in 
this article we saw that the non-antipodal version of the Sakai-Sugimoto model gives rise to 
the attractive potential as well as repulsive. Let us summarize our main results: 

• We argue that nuclear physics in the Nc ^ oc limit could be quite different from 
the real-life case of = 3, which limits the applicability of holography to zero- 
temperature nuclear physics. In particular, the ratio of kinetic to potential energy of 
the bulk nuclear matter scales with Nc as and consequently nuclear matter is 
a Fermi liquid for small Nc but becomes a crystalline solid for large Nc- We estimate 
the transition between liquid and solid phases happens for A^^c ~ 8, but this estimate 
is rather crude and should be taken with a large grain of salt. 

• Holographically, the attractive forces between nucleons arise from the coupling of the 
gauge fields living on the fiavor branes to the scalar fields parametrizing fluctuations 
of the those branes' geometry; in 5D, this coupling has form 



The antipodal Sakai-Sugimoto model has an accidental symmetry which 

forbids this coupling, so in that model there are no attractive forces. But the non- 
antipodal models don't have this symmetry, and consequently they do have the scalar- 
vector coupling ( 17.ip and hence the attractive nuclear forces. 

• At intermediate distances r between the two nucleons — larger than the nucleon radius 
p ~ Rkk/V^ but smaller than the Kaluza-Klein scale Rkk, see diagram on 
page [23] — both the attractive and the repulsive potentials have the 5D Coulomb 



eventually breaks down in the far zone r > {p/ y/T) ~ pV^ ~ Mj 




(7.1) 
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form, V{r) oc But the attractive potential has a smaller coefficient, so the net 

central potential is repulsive. 




1-C 



9 



-3 




(7.2) 



• There are similar attractive forces between different parts of the same baryon, and they 



• At longer distances r > Rkk ~ l/^mcsom the nuclear forces are dominated by 4D 
Yukawa forces due to the lightest meson with appropriate quantum numbers: The vec- 
tor isosinglet for the repulsive central forces and the scalar isosinglet for the attractive 
central forces. In the Sakai-Sugimoto model, the lightest scalar meson is heavier than 
the lightest vector meson, and consequently the attractive force has a shorter range 
than the repulsive force. It is not clear whether this un-realistic behavior is peculiar to 
the Sakai-Sugimoto models or a general problem of the large Nc limit. Indeed, in real 
life the lightest scalar meson is o"(600) but it's QCD origin is not clear. If it happens 
to be a bound state of two pions rather than a true qq meson, then in the large Nc 
limit this bound state will fall apart and the lightest surviving scalar meson would be 
heavier than the lightest vector meson a; (787). 

• The vector fields living on the fiavor branes are governed by the DBI + CS action, but 
usually the DBI part of the action is truncated to the lowest-order Yang-Mills terms. 



In the Sakai-Sugimoto model with A ^ 1, this approximation is valid at intermediate 
and long distances from baryons, but inside the instanton core of a baryon the non- 
abelian gauge fields become too strong to neglect the higher-order terms such as tr(J^^). 

We argue that the self-duality of the non-abclian magnetic fields saves the day and leads 
to approximate cancellation of the higher-order terms. This was known for instantons 
in flat space, and we show that this is also true for the Sakai-Sugimoto baryons; all 
we need to do is to slightly modify the self-duality condition for the non-abelian flelds 
to account for their coupling to the abelian electric and scalar flelds. The effect of this 
modiflcation on the baryon's radial proflle is quite small: 0(1/A) near the baryon's 
center and even smaller for r > p. 



reduce the baryon's radius by a factor (1 — Ca/^)^/^. Consequently, the isovector spin- 
spin and tensor forces between two baryons are reduced by an overall factor (1 — Ca/r)- 



>Cdbi oc str •>/det(^T27ra^ ^ •>/- det(g') x [Nf + 



^tr(^..^-)) 



(7.3) 
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We have also computed numerically the deviation of the baryon's profile from self- 
duality due to curvature of the flavor branes. The net deviation from a simple YM 
instanton is plotted on figure f l6.58p on page EH 

Our work gives rise to several open questions, the biggest of which is "What happens to 
the cr(600) meson in the large Nc limit?". The best answer for this question would be a 
lattice calculation of for several values of Nc, although it's not clear if such a calculation 
is possible with-present-day lattice sizes. (But thanks to Moore's law, it should be possible 
in a few years.) Alternatively, we can try several different models of holographic QCD and 
compare their predictions for the lightest scalar to lightest vector mass ratio. As of this 
writing, all known models have ratios > 1, with one exception [29] — but in that model, 
the lightest scalar meson is is a pseudo-Goldstone boson and probably does not couple to 
the other particles like the real a (600) meson. If future holographic models show the same 
pattern — the lightest true scalar meson is either heavier than the lightest vector or else is a 
pseudo-Goldstone boson whose couplings are suppressed — then most likely, in the large N^. 
limit of QCD there is no sigma meson and the attractive nuclear force has a shorter range 
than the repulsive force. But if we see a wide variation of the lightest-scalar-to-lightest- 
vector mass ratios between different models, then we wouldn't know what really happens in 
large- iVc QCD, but on the other hand, a holographic model with m^j < m^^ might also have 
a semi-realistic nuclear potential — repulsive at shorter distances but attractive at longer 
distances. 

We expect different holographic models to have different attractive / repulsive force 
rations at intermediate and short distances, and it would be interesting to see if any model 
has Ca/r- > 1- In such a model, net attraction between different parts of the same baryon 
would make it collapse to a singular point. Or rather, a classical baryon would collapse 
to a point, but quantum corrections would keep its size finite, perhaps 0{Rkk/ but 
much smaller than in the Sakai-Sugimoto model. Consequently, the net force between two 
such baryons would be repulsive at very short distances r < p Rkk/^ but attractive at 
intermediate distances p <^r <^ Rkk- At longer distances r > Rrk, the net force could be 
either attractive or repulsive, depending on the meson spectrum. 

Another open question concerns the dependence of the effective 5D action on the 't Hooft 
coupling A in the effective 5D action. In the Sakai-Sugimoto model, the flavor gauge 
coupling^ oc 1/A, but this power of A could be different in other holographic models. It 
would be interesting to find models where the flavor physics does not depend on A at all and 
to explore the baryons and the nuclear forces in such models. In particular, we would like 
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to see if such models have largish baryon radii, p ~ 1/Mmeson like in real life, rather than 
p ^ l/Mjncson as in the Sakai-Sugimoto model. For such large-radius baryons there would 
be no intermediate zone of distances; instead, the near zone (where two baryons overlap) 
would connect directly to the far zone dominated by 4D Yukawa forces. Consequently, for 
r ~ p both the baryon overlap and the curvature of the fifth dimension would be impor- 
tant, and the nuclear forces in this regime would be quite different from anything in the 
Sakai-Sugimoto model. 

There are also more general issues concerning meson spectra in holographic models. 
Apart from the specific mass ratios that are probably model dependent, there are general 
differences from the real-life mesons found in the Particle Data Book. For example, the 
5D scalar fields give rise to both scalar and pseudoscalar mesons in 4D (depending on the 
mode number in the fifth dimension), and their charge conjugation signs follow from parity, 
C-positive scalars and C-negative pseudoscalars. But in real life, all pseudoscalar mesons 
are C-positive rather than C-negative. 

Also, the high-spin mesons in holographic models have different physical origin from the 
low-spin mesons and consequently much larger masses. While the J = and J = 1 4D 
mesons are modes of 5D scalar or vector fields, the J > 2 mesons are semi-classical rotating 
open strings which start and end on flavor branes but at different points in space [33] . But 
in real life, both low-spin and high-spin mesons belong to the same Regge trajectories 

= a' X J + const (7.4) 

and there are no essential differences between them. 

We don't know what makes the holographic meson spectra so different from the real life, 
but it's almost certainly not the large Nc limit. It would be interesting to see if this problem 
is common to all holographic models — perhaps because it's inherent in the A — )■ oo limit — 
or if there are some model with more realistic meson spectra. If we can find such a model, 
maybe it would also have a light scalar meson and hence attractive net nuclear force at 
longer distances. 

Yet another open question concerns nuclear forces stemming from double-meson ex- 
changes, especially the double-pion exchange which produces a long-range attractive force. 
In holography, the double meson exchanges happen at the one-string-loop level while single 
meson exchanges happen at the tree level. This makes the double exchanges smaller by a 
factor 1/A, and also much harder to calculate. But if such calculation is feasible for some 



55 



holographic model, it would be very interesting to compare its result to the real-life nuclear 
force due to double-pion exchanges. 

Finally, there is a long-standing open problem concerning sensitivity of nuclear forces 
— and hence of the nuclear binding energy — to the pion's mass. Hopefully, holography 
can shed some new light on this old problem. Although holographic models usually have 
coincident flavor branes and hence zero current quark masses and massless pions, there are 
ways [361 133 EH |39] to explicitly break the chiral symmetry and give the pions a small mass. 
It would be interesting to see if a small but non-zero would affect the isoscalar central 
force between two nucleons, and whether such effect would happen at the tree level of string 
theory or only at the loop levels. 
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A Symmetrized Trace of the Non-Abelian DBI Action 

According to Tseytlin [M], the non-abelian version of the Dirac-Born-Infeld Lagrangian 
^J det{gMN + '^t^oi'^mn) works like this: First we focus on the spacetime indices of the 
Qmn + '^'^ol' ^MN and formally calculate the determinant of this dxd matrix while completely 
ignoring the gauge indices of the J^mn fields or the fact that they don't commute with each 
other. In other words, at this stage of the calculation, we treat each component Tmn as if 
it was a just real number rather than a generator of some non-abelian group. Second, we 
expand the square root of the determinant into a power series in the J^mn fields; again, we 
ignore the fields' non-commutativity and treat them as real numbers. Third, for each term in 
the expansion, we restore the gauge indices of the fields (in the fundamental representation 
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of a U{N) group), symmetrize the product of non-commuting fields, and take the trace 
str (j^M^NiJ^MiNi ■ ■ ■ J^AhNk^ = ^ tr ^all permutations of J^MiNiJ^M2N2 ■■■^i 



(A.l) 

Finally, we try to re-sum the power series in the J-" fields; if we a lucky, it might have a nice 
analytic form. 

In section 6, we had spherically-symmetric (in 4D) SU{2) fields ( I6.15p . and we had 
calculated the DBI determinant for those fields as 

det (g^n + 2Tia'J^mr) = det{g^n) x (l + aV^) (l + /^V^) . (iJHD 

In this Appendix, we calculate the symmetrized trace of the square root of this determinant 
and show that 

str V 1 + «2^) 1 + /32f2) = ■ fl6^ 

Clearly, expanding the square root on the LHS into powers of a and (3 produces all powers 
of = TiTi, so our first step is to symmetrize the product (r^)" with respect to all distinct 
permutations of the 2n Pauli matrices. 

Lemma: 

^ (2n-l)!! 

[(^^)Tsymm = (2^-1)!! distluct pcrmutatious of T^, T^, T.^ ■ ■ ■ ri„ri„ 

= (2n+l)xl2x2 (A.2) 

where 

(2n-l)!! = (A.3) 

is the number of distinct permutations of 2n matrices that come in n identical pairs. For 
the purpose of symmetrization, Xj and tj carrying different isovector indices i and j count as 
distinct; the summation over = 1, 2, 3 is done after the symmetrization. 

The lemma flA.2p is trivially true for n = and n = 1; indeed, for n = 1 there is nothing 



to symmetrize and [TjTj]_^ = TjTj = 6^ x 1 = 3 x 1. For the first non-trivial case n = 2, 



57 



the lemma works according to 



(r ) = TiTiTjTj 



symm 



^((TiTi)2 + {Ti,Tj} X TiTj^ 

^ ((3)2 + 25ij XTiTj = 9 + 2x3 
5 {i. e., 5 X 1). 



(A.4) 
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For higher n > 2 the simplest proof of the lemma is recursive. Let's group the {2n — 1)!! 
distinct permutations of the 2n matrices TiTiTjTj • ■ ■ into two sets according to the two left- 
most matrices being similar or distinct: {2n — 3)!! of the permutations start with TiTi (for 
the same i) while the remaining {2n — 2) x {2n — 3)!! permutations start with TjTj — or 
equivalently TjTi — with distinct i and j. Consequently, 



(2n-l)!!x[(fT] 



symm 



all distinct permutations of TiTiTjTjTkTk ■ ■ ■ 

(2n-3)!! 

Ti^i X permutations of TjTjTkTk ■ ■ ■ 

(2n-2)(2n-3)!! 

+ permutations of rjTjrfcTjkr^r^ • 

3x(2n-3)!! [(f')""'! 

L J symm 

+ Sij X (2n - 2)(2n - 3)!! \TiTj{T'^y^^ 
3(2n-3)!!x [(f^)"-'^ 



symm 



symm 



+ (2n-2)(2n-3)!! x [(r^y" i 
= (3 + 2n-2) X (2n-3)!! [(f^)""' 



symm 



symm 



and hence 



[(r-)»] 



symm 



X 



J symm 



(3 + 2n-2) X (2n-3)!! 2n + 1 



(2n- 1)!! 



2n-l 



(A.5) 



(A.6) 
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Applying this formula recursively, we get 



2n + l 



symm 2n - 1 



X 



symm 



2n + 1 2n-l 5 

X X ■ X 

2n - 1 2n-3 3 



2n + l 2n-l 

X X 

2n-l 2n-3 



symm 



(A.7) 



symm 



2n + l 



X 3 X 1 = (2r2 + 1) X 1. 



Quod erat demonstrandum. 

In terms of symmetrized traces, the Lemma tells us that 



str[(r-)T . tr([(r-)"] 



1 symm 



2(2n + 1), 



(A. 



which leads us to the following 

Theorem: the symmetrized trace of any analytic function 



(A.9) 



n=0 



of can be evaluated as 

[f{r')] =Y.CnX (str [{f'Y] = 2(2n + 1)) = Ux^ + 2 ) f{x) 



(A.IO) 



x=l 



In particular, the square root of the DBI determinant (I6.38p has symmetrized trace 
strv/(l + «27^)(l + /32f^) = (4«'^ + 4/3' + 2^ ^(1 + + ^2) 



2/Q2 



2 + 4a' + 4/3^ + 6a' 13 



(A.ll) 



as promised in eq. f l6.40p . 

To conclude this appendix, we note that the symmetrized trace (1A.11I) is bounded from 
above by the low-order tension+Yang-Mills limit 2 + 3(a' + /3^) and from below by its value 
2 + 6af3 for the self-dual fields, thus 



2 + 4^2 + 4/32 + 6^^2/32 
Va,/3>0, 2 + 6a/3 < , ^ h ^ h ^ 2 + 3 ^2 + 

v/(l + «2)(l + /3)2 - 



(A.12) 
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In terms of eq. fl6.4ip . these bounds amount to limits on P{a, (3), 
str v/(l + a2r2)(l + /32f2) = 2 + 6a(3 + P{a, /3)x{a-(3f, 0<P(a,/3)<3. (A.13) 
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